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ABSTRACT 



Large, fresh-water ice blocks (0.5m X 1.25m X 5m) were 
towed in Monterey Bay at speeds in range from 0.7 to 1.2 knots 
The objectives of the experiments included measurement of 
gross regression rates of ice surfaces, wake temperature, 
turbulent thermal boundary layer, ice interior temperature 
profile, and observation of shape changes over the melting 
period. The research was conducted over a period from 
October 1977 to January 1979. The measured regression rates 
at several points on the' ice blocks were compared to theoret- 
ical predictions using a turbulent flow ice ablation model 
developed by Dr. Owen Griffin of the Naval Research Laboratory 
Griffin's model predicted a regression rate of 280mm/hr 
compared with the measured value of 260mm/hr at a selected 
point. Temperature profiles in the ice interior and ice-water 
boundary layer -were used in a one-dimensional energy balance 
across the melting interface. Ice ripples were observed on 
all ice block tests. The ripples are important for both heat 
transfer and drag considerations. The wavelengths of the 
ripples were measured and the average wavelength determined. 
The measured average wavelength was 6.1 inches compared to 
the theoretical average predicted by Tatinclaux and Kennedy 
of 5.8 inches for a selected experiment. 
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INTRODUCTION 



I . 

A. JUSTIFICATION FOR STUDY 

There has been significant interest in recent years 
concerning the movement and deterioration of the large ice 
bodies that drift from the regions of the Arctic and 
Antarctic to the warmer waters of the Northern and Southern 
Hemisphere. The United States Navy in particular has been 
concerned about the impact these icebergs have on the 
temperature, density, and salinity of surrounding waters. 
Changes in these characteristics result in significant 
modification of the speed of sound through the water, which 
in turn complicates the problem of submarine detection, 
submarine navigation, and bottom navigation. The Navy and 
Coast Guard want the capability to predict with some degree 
of certainty the "life" of an iceberg and thus how it will 
affect surface ship operations over a period of time. In 
addition, as the need for additional sources of crude oil 
expands, offshore drilling in near-arctic regions surely 
will necessitate the movement of an iceberg by artificial 
means to avoid drifting into a drilling rig. How far to 
push or pull the iceberg, the power required, and the 
effect on the iceberg as a result of this movement are 
questions that must be dealt with before such an undertaking. 
A most important interest in iceberg deterioration, is 
the possibility of supplying water from icebergs to the 



17 



extremely arid regions of the world. A discussion of this 
area is given in Appendix A of this thesis. 
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II. BACKGROUND OF STUDY 



This thesis concerns a study conducted by the Naval 
Postgraduate School (NPS) in conjunction with the Naval 
Research Laboratory (NRL) on the ablation of ice bodies 
moving thorugh seawater. Relatively large (16' x4' xl.5') 
ice models were manufactured and moved through a body of 
salt water (Monterey Bay, California) . In addition, smaller 
models were placed in a laboratory flow channel for experi- 
ments and analysis. The ultimate goal was the ability to 
predict the overall ablation rate of an ice body moving 
through seawater. The prediction should be verified by 
comparing the experimental results for this thesis with 
ablation theory previously published by Owen Griffin [1,2] 
of NRL. 

A first important step in the study was to analyze the 
scaling problem relating a large iceberg (15,000' x2500' x700') 
to a much smaller model. For the tests in Monterey Bay, 
the size of the model was constrained by manufacturing, 
crane, and transportation facilities available. In addition, 
for the experiments in the flow channel, the model size 
was limited by the flow channel equipment at NPS. The 
scaling is important because the experimental results are 
more valuable if the information can be extrapolated to a 
full-scale iceberg. The physical effects of the ice moving 
relative to the seawater (i.e.: shape change, ripple 
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effect, and wake temperatures) were examined. The thermal 
boundary layer between the ice and water was measured and 
compared with that theoretically predicted. Finally the 
heat transfer through the ice was analyzed by comparing 
temperature profiles in the ice over a period of time with 
those predicted from a solution of the equations governing 
transient heat transfer. 
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III. ICE MODEL SCALING AND INITIAL LABORATORY TESTS 



A . GENERAL 

The problems encountered in modeling a large ship for 
laboratory water channel experiments are also inherent, 
but to a greater magnitude, in modeling a large ice body. 
This is due to several factors, the most predominant being 
the immense size of the ice bodies to be modeled coupled 
with the low speed of transport (1.0 kt) . Also resolving 
the hydrodynamic requirements with those of the heat 
transfer analysis to realize an acceptable model from which 
melt rates can be extrapolated to the full size ice body is 
also an important consideration. Water channel and towing 
tank modeling in general are discussed thoroughly in Chapter 
7 of ref [3], and modeling particular to ice bodies is 
discussed by Murphy [4] and Griffin [1,2]. Murphy dis- 
cussed the methods used to determine towing forces and mass 
loss during transit, and Griffin treated the heat transfer 
problem specifically to predict ablation rates of large 
ice bodies. 

B. HYDRODYNAMIC REQUIREMENTS 

Two well-understood dimensionless groups of variables 
are commonly used to relate the forces on a body moving 
through a fluid. These are: 



Reynolds number (Re) 



VL _ inertial effects 
v viscous effects 
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Froude number (Fn) 



V 



inertial effects 
gravity effects 



( 2 ) 



/gL 



The Reynolds number is important in defining whether a 
fluid over a body is laminar or turbulent, a fact that has 
implications in both towing force and heat transfer analyses. 
Since the kinematic viscosity (v) is a property of the 
fluid, the only variables for Reynolds number which can be 
easily controlled by the experimenter in ice model testing 
are the velocity of the fluid (v) and the characteristic 
length of the model (L) . For ideal modeling, the Reynolds 
numbers of the model and prototype should be identical to 
assure the same laminar or turbulent flow conditions. This 
means that a much smaller model must move through the fluid 
at a very high velocity. The high speed of the test fluid, 
however, is in direct conflict with the Froude number which 
is used in determining residuary resistance of a prototype. 

In addition, the Froude numbers of the model and prototype 
should be identical for ideal modeling. Note that the 
Froude number also has only fluid velocity and model length 
as independent parameters. 

The solution commonly used in ship modeling for the 
conflicting Reynolds number - Froude number requirements 
is to scale the velocity and model size to satisfy Froude 
number similitude. Reynolds number scaling is unavoidably 
violated. However, models are tested at either moderately 

g 

high Reynolds numbers (10 ) to assure natural transition to 
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turbulent flow over forward regions of the model or with 
the use of artificial boundary layer trip mechanisms placed 
so as to stimulate transition at the predicted natural 
transition point on the model. 

C. HEAT TRANSFER ANALYSIS REQUIREMENTS 

Griffin [1,2] at NRL proposed a heat transfer model 
for a large glacial ice body moving through seawater with 
turbulent flow. In Griffin's model the local turbulent 
melting rate is given by: 




DSt 



h , TURB 



(3) 



where : 



^TURB 



Local turbulent melting rate (m/s) 



D 



Stefan number parameter determined bv: 



D 



A 



(4) 



1 + B 



and: 



Seawater Stefan number = Cpp(T - Tp)M 



-1 



A 



Ice Stefan number - (T„ - T )M 

For O 



-l 



B 



St. 



Local turbulent Stanton number 



(5) 



h , TURB 
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where R is the ratio of interfacial heat transfer with 
and without melting; R is assumed to have the same value 
in laminar as well as turbulent flow. 

The local turbulent Stanton number on a flat plate 



is : 



St 



h , REF 



Nu 

x 

Re Pr 
x 



.0297 Re 



-1/5 



x 



1 + 1.48 Re 1//10 Pr 1//6 (Pr-1) 
x 



( 6 ) 



By substituting eq. (6) into eq. (5), and then into eq . 
(3), the variables affecting the local turbulent heat 
transfer rate can be seen as: 



.0297 Re 



-1/5 



^TURB DR [ 



1 + 1.48 Re ~ 1/10 Pr~ 1/6 (Pr-1) 



(7) 



Since Prandtl number (Pr) and Stefan number (D) are 
properties of the seawater and ice, both numbers are the 
same for model and prototype. Likewise R, which is depen- 
dent on the properties of ice and the solid-liquid phase 
transformation, is also invariant between model and 
prototype . 

For ideal modeling: 

^model _ ^ ^ g 

^prototype 



From eq. (7) it is obvious that agreement between Reynolds 
numbers of the model and prototype is paramount for ideal 
modeling as stated by eq . (8). 
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From the above consideration, model testing of ice 
bodies is limited to a regime of Reynolds numbers and 
Froude numbers as shown in Figure 1. The points indicated 
on the figure are for experiments that will be discussed 
later in this chapter and in succeeding chapters. 

D. WATER CHANNEL FACILITIES AT THE NAVAL POSTGRADUATE SCHOOL 

The Naval Postgraduate School, has several water channels 
which are pictured in Figures 2-5. The characteristics of 
each channel are summarized in Table I. Since buoyancy was 
considered an important factor in the heat transfer modeling 
analysis, the water tunnel shown in Figure 2 was eliminated 
as a testing facility. The water channel in Figure 3 was 
also eliminated because of the extensive modifications 
required to obtain a uniform velocity profile over the 
width of the channel. The remaining two water channels 
were utilized; however, their relatively small size pro- 
hibited the achievement of Reynolcs numbers and accompanying 
turbulent flow conditions equivalent to those on a large 
ice body. The inability to achieve these conditions was 
a direct result of the conflicting requirements of high 
water channel velocities (2-4 ft per sec) and long model 
lengths (3-6 ft) necessary for moderately high Reynolds 
numbers and the low velocities needed for satisfactory 
Froude number scaling. 

Large water channels, similar to the new Environmental 
Protection Agency water channel in North Carolina, with 
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dimensions of L = 80 ft, W = 8 ft, and D = 4 ft, would 
allow the testing, of models large enough to achieve a more 
satisfactory scaling of Reynolds number, geometric size, and 
velocity. This scaling would be achieved by virtue of the 
relatively long models (8-16 ft) that could be tested at 
the relatively low velocities (0.1 -0.5 ft per sec) required 
However, the time, expense, and coordination involved in 
conducting experiments remote from NPS were deemed prohibi- 
tive; therefore the decision was made to conduct a series 
of large scale model tests in Monterey Bay adjacent to 
NPS. The details of the experiments in the bay are dis- 
cussed in Chapters IV-VI . 

A review of literature dealing with large ice bodies 
revealed no experimental data on melt rate, shape change, 
or regression rates of fresh water ice moving through sea- 
water. Therefore, in order to gain at least a qualitative 
understanding of what these characteristics were, a series 
of model tests were conducted in the NPS water channels. 

E. INITIAL LABORATORY EXPERIMENTS 

The water channel selected for the initial tests is 
pictured in Figure 4 with characteristics listed in Table I. 
The size of the ice models was selected to avoid both 
blockage and side and bottom wall flow interference in 
the channel. 

Four rectangular pieces of ice were constructed with 
overall dimensions of L = 8.0 in, W = 1.5 in, and D = 1.3 in 
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TABLE I 



Naval Postgraduate School Water 
Channel Characteristics 



Figure No. 


Test Section 
Dimensions 
L, W, H (in) 


Velocity 

Range 

ft/sec 


Re 


2 


16/4/8 


0-30 


10 


3 


96/48/5 


LO 

O 

1 

o 


10 


4 


36/12/16 


0. 2 - 3.0 


10 


5 


36/24/8 


o 
1 — 1 

1 

o 


10 
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thus achieving a satisfactory ratio of channel width to 
model width of 10:1. The models were frozen in a standard 
freezing compartment at a temperature of 10 °F. A wire 
restraining device was frozen into the body of the model for 
handling during testing. To test the effect of a stream- 
lined shape during identical flow channel conditions, one 
model was constructed with all sides sloping at 30 degrees 
from the vertical. Overall dimensions at the top surface 
were L= 9.7 in, W = 3.7 in, and D = 1.2 in. The ratio of 

channel width to model width for this model was a less 

satisfactory 4:1; however, the low channel water flow 
velocity virtually precluded any significant wall inter- 
ference . 

Each model test was similar and consisted of transferring 
the model from the freezing compartment to a scale immedi- 
ately above the flow channel where the model was weighed 
and photographed simultaneously in the side view. The model 

was lowered inot the flowing water and allowed to float 

naturally. Directional control of the model was accomplished 
by attaching the model restraining device to a cross-beam 
over the channel. At regular intervals (1 min) , the model 
was removed, weighed, and photographed until the restraining 
device separated from the model. Dyes were frozen into the 
models, or, alternatively, were injected in front of the 
model for visualization of the flow patters. In four out 
of five tests conducted, the data on weights of the receding 
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models, air and water temperature, flow velocity, and 
regression rates at several areas of the models were noted. 

The total time in water for each model was as follows: 

(1) 2 min/50 sec, (2) 2 min/30 sec, (3) 3 min/30 sec, 

(4) invalid (due to hookup difficulties), and (5) 2 min/30 sec. 

F. WATER CHANNEL TEST RESULTS 

The first four models melted through a series of shapes 
similar to those illustrated in Figure 6 . The second model 
tested is shown in Figure 6; results of the melting are 
given in Table II. The rates of regression at three 
selected surfaces for the model in Figure 6 are given in 
Table III. 

The frames of Figure 6 are side views with the front 
of the ice oriented to the right of the photograph. The 
streamlined ice model was the fifth to be tested. The 
sequenced melting shape configuration for this model is 
shown in Figure 7 with similar orientation. It was noted 
that despite the significant geometric differences between 
the streamlined model and the four previous models, the 
melting configuration is practically identical. The regres- 
sion on the front, bottom, and rear of the models can be 
seen in Figures 6 and 7 which clearly illustrate the area 
of highest regression as that on the bottom, approximately 
three-tenths of the model length from the front. The 

laboratory conditions included relatively small Reynolds 

4 

numbers (approx. 10 ) compared with the large numbers 
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TABLE II 



Melt of Second Rectangular Ice Model 
(See Figure 6) 



Photo Frame 



Weight (g) 



Time (sec) 



Upper 

Upper Middle 
Lower Middle 
Lower 



152 


0 


92 


60 


49 


120 


15 


150 



Air Temperature - 23 °C 

Flow Velocity 



Water Temperature 21.3°C 
0.55 f t/sec 



TABLE III 

Regression Rates of Second Rectangular Ice Model 

(See Figure 6) 

Average Rates for 120 Seconds 



Rear 



Bottom (neck) 



Front 



0.10 



mm 

sec 



0.16 



mm 

sec 



0.13 



mm 

sec 



30 



9 

(approx. 10 ) associated with large Atlantic ice bodies. 

Given these small Reynolds numbers achieved in the labora- 
tory, it was concluded that the flow conditions in the 
laboratory were largely laminar and that the decisive necking 
effect observed was associated with the reattachment of the 
laminar boundary layer initially separated from the model 
by the sharp corners of the blunt bow. 

For the melting observable in Figures 6 and 7 along the 
profile surfaces of the models, the Nusselt numbers, which 
measure the effectiveness of heat transfer, can be calcu- 
lated to show this effectiveness at any distance along the 
profile surface. Figure 8 shows a plot of Nusselt number 
as a function of length along the ice model. The plot 
shows a well-defined region of intense heat transfer at a 
distance of approximately 2D (D = Depth) from the front of 
the model and an accompanying decay of melting at approximately 
3.5D from the same point. A representative theory for 
laminar flow on a flat plate with constant surface tempera- 
ture is also shown on the plot in contrast to highlight the 
empirical results observed in the laboratory. 

To explain the high melting rate in the forward third 
region on the models, the following hypothesis is proposed: 

Due to the initially sharp corner on the model, the flow 
separates from the model as it passes the corner point. 

In the region 0 - 2D, there occurs a recirculating flow with 
reattachment near X = 2D. As is well known, the area where 
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reattachment occurs is a region of high heat flux as analyzed 
by Meyer [5] and Koseman [6] . Figure 9 illustrates this 
hypothesis. 

An alternate hypothesis involves a trapped vortex 
filament as shown in Figure 10. The sense of the vorticity 
is shown by the arrows and follows the customary right-hand 
rule. The vorticity originates because of the flow separa- 
tion from the blunt bow. Initially the vorticity is shed; 
but as the ice ablates, the vortex filament becomes trapped. 
Due to the vorticity, the warm channel water is circulated 
close to the ice resulting in large temperature gradients 
and accompanying high heat flux. The observations made during 
the model experiments were consistent with either the 
reattachment or trapped vortex filament hypothesis. 

G. CONCLUSIONS 

The normal problem of resolving the conflicting require- 
ments of Reynolds number and Froude number scaling in 
modeling ships is magnified when applying standard modeling 
techniques to large ice bodies. This difficulty in modeling 
is primarily due to the extremely large size and low trans- 
port speed of these full-size bodies. The high Reynolds 
number required for the model is necessary to assure turbu- 
lent flow and to minimize the difference in theoretical 
melt rate between the small-scale ice model and full-size 
ice body. The only practical controlling variables available 
in developing the model are the size of the model and the 
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velocity of the water. The flow separation around the 
frontal corners of. small models results in an extremely 
high heat flux and accelerated melting rate at the point 
of reattachment. This severe melting tends to mask all 
other characteristics that could be determined from such 
a small model. The test results indicated that large scale 
model testing in Monterey Bay, despite the lesser controlla- 
bility of environmental conditons, was needed in order to 
determine more accurately the melt-rates, shape changes, 
and other characteristics of a full-scale ice body underway 
in seawater. 
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IV. LARGE-SCALE ICE MODEL EXPERIMENTS 
IN MONTEREY BAY 



The following is a discussion of the ice model size 
and shape determination, freezing container design, freezing 
technique, and general test procedure and a summary of 
each test conducted in the Bay. Specific data results 
and theoretical correlations from the various tests are 
presented in Chapter VI. 

A. SIZE AND SHAPE DETERMINATION OF MODEL 

Several major considerations of both scientific and 
practical nature were important in determining the minimum 
and maximum size of the ice body model to be tested in 
Monterey Bay. The controlling minimum dimension was the 
model length, which must be sufficient to achieve natural 
transition to turbulent flow at speeds of 1-2 ft/sec (Re 

g 

10 ) . The model required sufficient volume to withstand 
the following test conditions: a one-hour loading and 

transit time (20 mi) from the vendor to the launch site on 
an open truck with minimal insulation, a second hour period 
at the launch site for testing preparations and launch, a 
30-minute underway period in the Bay, and recovery from the 
water intact for examination of resulting changes in physical 
characteristics. The maximum size of the model was limited 
primarily by the size of the freezing compartment which 
was the only large compartment available in the immediate 
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geographical area. However, other factors such as the 
load capacity of a flatbed serai-truck, pier side space 
restrictions, the available crane lifting capability, and 
the maneuverability of a boat with ice in tow within the 
restricted harbor were also important in determining the 
maximum size. The utmost consideration throughout the 
development of this previously untried test series was the 
complete safety of all personnel during the tests. 

The model shape was desired to simulate a tabular 
Antarctic iceberg similar to those that continually break 
off from the ice shelves and drift Northward. These ice- 
bergs have considerable length as compared to width and 
depth and thus have good longitudinal hydrodynamic stability 
for towing. The ends of these icebergs are blunt and have 
a flat table-top appearance above the waterline. The 
tabular icebergs are readily duplicated in shape using an 
easily manufactured rectangular container. 

The result of the foregoing size and shape considera- 
tions was an ice model with dimension of L = 16 ft, W = 4 
ft and D = 1.5 ft. This particular size and shape satis- 
fied all of the requirements and was used unchanged through- 
out the entire series of experiments conducted in Monterey 
Bay . 

B. CONSTRUCTION OF THE ICE MODEL FREEZING CONTAINER 

The container was manufactured by a local vendor from 
galvanized steel sheeting with welded seams. The 
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specifications for the container are listed in Table IV, 
and the container is pictured in Figures 11 and 12. Suffi- 
cient integral bottom support structure was included to 
facilitate the lifting of the 6000 lb. ice model and test 
accessories. Removable plugs were installed to allow flooding 
of seawater into the container during launching and draining 
of this water after launch of the. model. This simple design 
proved successful requiring no modifications throughout the 
test series. 

C. ICE MODEL CONSTRUCTION 

For the first two experiments conducted in Monterey Bay, 
the container was placed in the vendor's freezing compartment, 
filled to the 18-inch level with demineralized water (18 ppm 
dissolved solids), and frozen as a single layer. Prior to 
filling with water, the towing bollard, seen in Figure 13, 
was suspended over the container centerline, 48 inches 
behind the front wall, with the bollard base 9 inches above 
the bottom for frrezing into the ice. The bollard position 
was determined to afford maximum lateral towing stability, 
while remaining intact in the ice during the course of the 
test with rapid melting of the front and bottom surfaces 
of the ice. This position also precluded any hydrodynamic 
interference with the water flowing past the ice and the 
subsequent masking of natural physical characteristics by 
this flow interference. 
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TABLE IV 



Specifications of Ice Model Container 

1. Dimensions (bottom) 

Length - 16 ft 
Width - 4 ft 

Height - 1.66 ft 

2. Slope of Sides 

Seven degrees out on all sides from bottom. 

3. Material 

Sides - 12 GA galvanized mild steel. 

Bottom support structure - mild steel of 
sufficient strength to support 6500 lbs. 

4. Drain Plugs 

Install two 2" dia drain plugs on center line 
bottom of front and rear sides. 
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The size of the bollard vertical post and the horizontal 
base plate were designed to minimize the towing stress 
concentrations on the ice. This was accomplished by dis- 
tributing the towing forces over a wide area in both 
horizontal and vertical planes, without significantly 
weakening the ice at the critical vertical lateral lane 
passing through the bollard installation point. The bollard 
protruded above the ice, as pictured in Figure 12, and also 
served as an anchoring point for instrumentation during the 
tests . 

Two problems developed from freezing the ice as a single 
18-inch layer. First, the expansion of the ice during 
freezing resulted in a large upheaval of the top surface 
as the core solidified. The upheaval of the ice was accom- 
panied by large vertical fissures up to one inch in width 
on the surface and extending several inches into the ice 
at some points as shown in Figure 12. Second, as a result 
of the upheaval, the effective freeboard (distance from top 
surface to waterline) was significantly lowered allowing 
seawater to flow over the top surface of the ice. The result 
was rapid erosion of the upper surface of the bow and 
fissured areas with premature disintegration of the model 
as shown in Figure 14. 

To eliminate the model's top surface upheaval and low 
freeboard, two modifications were made to the model construc- 
tion technique. First, the ice was frozen in a series of 
layers 1-2 inches in depth to minimize the extreme expansion 



38 



and subsequent fissures. Second, a hollow wooden box 
(Archimedes box) was frozen into the ice immediately behind 
the bollard and at the same nine-inch depth. The Archimedes 
box, which is identified here as the core box, reduced the 
weigh t-to-volume ratio of the model and increased its free- 
board 1.5 inches. The core box and the installation are 
shown in Figures 15 and 16. The combination of layered- 
freezing and flotation box installation resulted in a 
hydrodynamically stable, non-f issured, model with sufficient 
freeboard for testing as shown in Figure 17. 

The addition of a small bollard on the rear of the core 
box improved the waterborne handling of the model. Further- 
more the box served as a reference centerline and anchoring 
point for test instrumentation. The only further modifica- 
tion made to the model throughout the remainder of the test 
program was to include the towing bollard as an integral 
part of the core box as shown in Figure 18. The revised 
method of holding the bollard resulted in a weight savings 
of approximately 25 lbs, while reducing the quantity of 
foreign material frozen into that area of the model. From 
all considerations ranging from transportability to hydro- 
dynamic similitude, the model was satisfactory for the 
test series. 

D. GENERAL TEST PROCEDURE 

The test procedure described below is typical of the 
series of experiments. Specific deviations from this 
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general procedure are described in the summary of each 
experiment later in this chapter. 

Approximately two weeks prior to each experiment, the 
ice container was delivered by the Naval Postgraduate School 
Public Works Department to Growers Ice, Salinas, California. 
The container was placed in a freezing tunnel and filled 
with a 1 - 2 inch layer of demineralized water and allowed 
to freeze. Once the layer was frozen, the tunnel was opened 
and allowed to warm to outside ambient temperature, and an 
additional layer of water was added over the frozen ice 
for freezing. This process was continued to the 9-inch 
level at which time the core box, bollard, and thermocouples 
were installed. The process was then continued to the 
18-inch level with the model components firmly frozen in 
position. The average temperature of the tunnel during 
freezing was 0°F, and the total time to complete the 
freezing process was 10 -12 days. The model was further 
cooled to minus 18 - 20 °F twenty-four hours prior to pickup 
in order to inhibit melting during the transit and testing 
preparation period as previously discussed. 

In order to exploit the early morning quiesent wind and 
sea conditions, to avoid marine traffic interference, and to 
minimize the melting by direct sunlight, the tests were 
conducted between 0700 and 1100. The ice was brought to the 
launching site at the U.S. Coast Guard Pier, and launch 
preparations were made. Launch preparations included 
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inspection of ice for fissures, dimensional measurements, 
hookup of instrumentation, and attachment of towing cables 
and lifting rigging as shown in Figures 11, 12, 16, and 18. 

At launch time the entire assembly was weighed and lowered 
into the water as pictured in Figures 19-22. As seawater 
flowed into the container, the assembly was lowered slowly 
by the crane until the buoyant ice model floated freely as 
seen in Figure 23. Connections were made to the towing boat, 
the model pulled clear of the rigging, and the underway 
portion of the experiment commenced as pictured in Figure 24. 

Primarily for safety purposes, considering that the 
tests were to be conducted in a restricted harobr, the 
towing vessel for the first experiment was a Coast Guard 
patrol vessel as shown in Figure 25. All subsequent tows 
were conducted utilizing a Boston Whaler with a small out- 
board motor pictured in Figure 26. A wire bridle connected 
both arms of the towing bollard to a tow rope, and the 
thermocouple wire bundle was streamed parallel to the tow 
rope as illustrated in Figure 27. The ice was towed at an 
average speed of 1 knot (1.78 ft/sec); the speed was 
determined using a Marine Advisor Ducted Current Meter and 
Readout, shown in Figure 28, which was suspended on a rigid 
arm from the towing boat. The model was tewed for approxi- 
mately thirty minutes in a "race track" shaped pattern 
approximately 600 yards (.34 mi) in length on a side. The 
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actual time of each tow varied, with the object being to 
have the model underway for a sufficient period to allow 
the full development of surface ablation characteristics, 
while maintaining adequate structural strength for recovery 
at the pier. During the tow, temperatures of the ice 
interior, ice-seawater boundary layer, and model wake were 
recorded. 

Recovery was accomplished by positioning the ice model 
in a web cargo net with a tubular-steel bottom frame and 
raising the entire assembly onto the pier as sequenced in 
Figures 29-31. The assembly was weighed as before to 
determine the total mass lost during the test. Once on the 
pier, the model was photographed thoroughly, surface ablation 
characteristics noted, and overall dimensions recorded. 
Subsequently the ice was broken for disposal and recovery 
of the core box and thermocouples. 

E. SUMMARIES -OF TESTS IN MONTEREY BAY 

The following sections are brief discussions of each 
test conducted in the Bay. The discussion includes details 
of deviations from the general test procedure and evolution 
of the test procedure as experience was gained on each 
experiment. 

1 . Test No. 1, September 23, 1977 

The ice was transported from the vendor to the launch 
site after a four-day freezing period as a single large 
block. The characteristic "hump" and fissures previously 
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described were present as a result of the freezing technique. 

A large wooden reference grid was attached prior to launch 
as seen in Figure 12; the grid provided a reference for 
measuring regression of the ice surface. The container was 
placed in the water, and the ice floated free and was taken 
under tow as pictured in Figure 32. The large towing vessel 
had difficulty maintaining the desired speed of 1 knot. 

The low freeboard and numerous fissures allowed water to 
flow over the top surface resulting in premature breakup of 
the ice and detachment of the towing bollard while underway. 

The bollard and reference grid were recovered; 
however, the ice was not recovered. Therefore, significant 
data were not obtained. Though this test was largely incon- 
clusive from a scientific view, the photographs obtained and 
the experience gained in the launching and towing procedure 
were useful for future experiments. 

2 . Test No. 2, September 26, 1977 

The ice was removed from the freezing compartment 
after a three-day freezing period as a single layer as before. 
However, three days were apparently insufficient to solidify 
the model completely. Ice formed on the outside with water 
on the inside. The water was covered by a thin layer of ice. 
The thin layer ruptured prior to arrival at the launch site 
as illustrated in Figure 33. The overriding point which is 
apparent from Figure 33 is that the center of the ice did 
not have sufficient time to solidify. The second experiment 
was scheduled only a few days after the first experiment in 
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order to obtain additional information for presentation at 
the First International Iceberg Utilization Conference on 
October 3, 1977. 

The container with the partially melted ice was 
returned immediately to the vendor to resume freezing. 

3 . Test No. 3, October 3, 1977 

The ice model for Test No. 3 was frozen during a 
five-day period as a single layer. Note that the ice for 
Test No. 3 was in part the same as for Test No. 2. The 
"hump" and fissures were again present on the top surface. 

The reference grid was attached, and the model was launched. 
The model freeboard was slightly higher as pictured in 
Figure 34 but was still insufficient to prevent water flow 
over the top surface. The towing vessel used was a 16-foot 
Boston Whaler with an outboard motor which had little diffi- 
culty maintaining the desired speed. Underway photographs 
were obtained; however, the reference grid could not be 
used to determine regression rates over the length of the 
model because of the relatively small differential in these 
rates and the inability to obtain photographs directly over 
the ice . 

The ice model was recovered successfully for the 
first time as illustrated in Figures 35 and 36. The recovery 
net shown in Figures 35 and 36 did not have the bottom 
tubular steel frame pictured in Figure 30. The tubular 
frame, which was included for later tests, spread the net 
over the width of the model. Also, the frame submerged 
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the net in the water to ease the positioning of the model 
in the net for recovery. In addition, the frame distributed 
the lifting forces over the bottom of the model, thus 
minimizing loads applied to the delicate sides and corners. 

The model was placed on the pier where the weight, 
dimensions, and surface ablation characteristics were 
recorded. Numerous photographs were taken similar to 
Figure 37 which illustrate the erosion caused by the sea- 
water flowing through the fissures on the model. Two signi- 
ficant surface features were noted: First, a heavily eroded 

cavern-shaped area on the bow as shown in Figure 38 and, 
second, as illustrated in Figures 39 and 40, a series of 
well-defined, vertical, periodic ripples extending the length 
of either side of the model. These two surface character- 
istics also were observed in subsequent experiments and are 
discussed with other test results in Chapter VI. 

One surface characteristic that was not present on 
this and subsequent tests was the decisive "necking" observed 
on the laboratory models. However, if this process did 
occur, it was only present in the first few inches and then 
was so insigifniciant that its effect was masked by the 
ripple system on the model. An amplifying discussion of 
this surface characteristic is presented in Chapter VI. 

This test was judged successful in that a complete cycle 
from launch to recovery was accomplished along with obser- 
vations of the ice model dimensions and weight which were 
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measured as a result of the cycle. However, modifications 
to the freezing technique, increased freeboard, and improve- 
ment of the useless reference grid were necessary for the 
next model test. In addition, temperature sensing devices 
and ripple measuring techniques were required to increase 
the amount of experimental information obtained. 

4 . Test No. 4, November 6, 1977 

Test No. 4 marked the first attempt to obtain com- 
prehensive data on the ice interior temperatures, ice- 
seawater boundary layer temperatures, and ripple wavelength. 
The model was frozen using the layered technique previously 
described, with the bollard, core box, and thermocouples 
frozen at the nine-inch level. Figure 41a is a top view of 
the model illustrating the placement of these devices with 
an amplifying pictorial view in Figure 42. The numerals in 
Figure 41a denote the amplifier channel for that particular 
thermocouple . 

The thermocouples were 1/16 inch diameter Omega 
fast responses Copper-Constantan Type T with beaded sensory 
tips. Five of these thermocouples were bonded together near 
the tips to form an array of five thermocouples spaced at 
1/16 inch intervals (channels 11-15) . The array was mounted 
on the bollard (Figure 41b) and placed in the ice with the 
outermost thermocouple 1.5 inches from the side surface. 

The purpose of the array was to indicate both the interior 
ice temperature profile as the side receded towards the 
array and the ice-seawater boundary layer profile as the 
array became exposed to seawater. 



Interior ice temperatures were recorded on channels 
16-20 (Figure 41a) , The output of a particular thermcouple 
was selected by a channel selection and amplified with 
accompanying electronic reference junction (zero degrees 
Centigrade) as shown in Figure 41b. The output of the digital 
voltmeter was recorded by voice using a tape recorder. The 
data were later transcribed and converted to °c using the 
thermocouple conversion tables. 

The model was launched without incident and, as 
pictured in Figures 43 and 44, had a completely flat top 
surface and sufficient freeboard. At regular intervals all 
thermocouples were read and recorded. The model was towed 
until all five thermocouples in the array indicated ambient 
seawater temperature and significant undercutting was observed 
as in Figure 45. At that point the model was recovered 
without difficulty. Once clear of the water, the extent of 
the differential in melt rates between submerged and above- 
water surfaces readily became apparent as shown in Figures 46 
and 47. Close examination of Figure 47 shows the exposed 
array protruding down from the top surface. 

Once on the pier, the various layers comprising the 
model were very well defined as shown in Figures 48 and 49, 
as were the ripples and heavily eroded bow section noted in 
the previous experiment. The ripples were measured by 
pushing a multi-fingered contour gage against the surface 
until the fingers conformed to the surface shape. The 
contour was traced onto a strip-chart paper in succession 
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with the end result being a profile of the ice surface 
drawn on the paper. Relative depths were measured using a 
metal reference frame prepositioned on the centerline of 
the core box as illustrated in Figure 50. 

Test No. 4 was considered successful, although the 
ripple measurement method was cumbersome and needed improve- 
ment. The effect of the ice layering on the regression rates 
and ripple formation was a question. In the next test 
the freezing layer thickness was increased to three inches 
to evaluate any dependence on ice layer thickness. In 
addition, the rapid regression at the ice-seawater boundary 
layer with the subsequent exposure of the array resulted in 
difficulty of monitoring these channels with the current 
recording equipment. For future tests a single thermocouple 
with direct output to a portable strip-chart recorder was 
utilized to collect boundary layer temperature profile 
data. 

5 . Test No. 5, December 10, 1977 

The alternative 3-inch layer freezing technique was 
applied to this model with adverse results . The increased 
net expansion of the ice buckled the bottom and the sides 
of the core box breaking water tight seals . Broken seals 
allowed water to flow into the box with subsequent freezing. 
As a result the buoyancy of the box was reduced. The top 
of the core box was opened at the launch site, and the ice 
within the box was chipped out as shown in Figure 51. The 
box was resealed, and the model was launched. Although the 
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bottom surface of the box was buckled upward reducing volume, 
sufficient buoyancy remained to provide an adequate free- 
board during the test. 

The thermocouple arrangement for the ice model is 
shown in Figure 52. The arrangement was designed to measure 
boundary layer profiles near the bow stagnation point (Channel 
13) and at the indicated side points (Channels 11 and 12) . 
Since only one portable strip-chart recorder was available, 
the thermcouples were placed accurately at carefully selected 
distances from their respective receding edges so that they 
would be sequentially exposed and recorded. The sequence 
of channels was 13, 11 and 12. The sequence was determined 
using previous regression rate data. The thermocouples were 
monitored initially on a digital readout and then trans- 
ferred to the strip chart recorder when the output approached 
0°C. Each channel was switched back to a digital readout 
when the strip chart indicated ambient seawater temperature. 

In addition to the boundary layer, an effort was 
made to measure the temperature profile across the wake of 
the model. The purpose of this effort was to determine the 
extent of cooling experienced by the water flowing past the 
ice. If cooling proved to be significant, a second ice 
body towed in the wake of a first would be "insulated" 
effectively by cooler water; and the melting rate of the 
second body could be reduced considerably. The temperature 
measurement was accomplished with an array of seven thermo- 
couples positioned across a bridge between two pontoons as 
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shown in Figures 53 and 54. The pontoon bridge was pulled 
behind the ice by a rope attached to the rear bollard on 
the model. Figures 55 and 56 illustrate the towing arrange- 
ment which had the thermocouples positioned one model width 
(4 ft) behind the ice. A distance of 4 feet was a compromise 
based on being sufficiently close to the model to detect 
cooling while remaining clear of the immediate eddies behind 
the model which could distort overall temperature gradients 
in the wake. The array measured the wake temperature at 
a depth of 6 inches, and over a lateral distance 30 inches 
on either side of the centerline. Thermocouple output was 
voice recorded from the digital indicator at periodic 
intervals during the test. 

During the underway portion of the test, the added 
drag of both the pontoon bridge and fouling by floating 
kelp proved beyond the capacity of the electric outboard 
motor to maintain desired speed. Thus only 80 percent of 
desired speed was attained during the test. The lower speed 
had the effect of decreasing the ablation rate and increasing 
the time required for the boundary layer thermocouples to 
become exposed. Although the boundary layer thermocouples 
eventually became exposed and their outputs recorded, 
recovery of the model was necessary due to a large crack 
forward of the core box precluding exposure of the stagnation 
point thermocouple. 

The recovered ice is pictured in Figure 57 where 
close examination shows the stagnation point thermcouple 
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tip visible through the ice surface, while the forward side 
thermocouple tip clearly protrudes thorugh the side of the 
ice. The crack forward of the core box also can be seen 
clearly in Figure 57. The bow erosion, ripples, and under- 
cutting were again observed in this test and are pictured 
in Figures 57-59. These features were less pronounced, 
however, due to the reduced towing speed. Ripple data were 
obtained using the same contour gage method. An improved 
method was under development but not available for Test No. 

5. The thicker freezing layers also were well defined, but 
their effect on the ripple system and undercutting was 
inconclusive due to the lower speed of Test No. 5. 

A matter of practical importance became apparent 
during this test in that the complexity limitations of the 
ice model testing were being rapidly reached. The towing 
vessel which required an operator, line handler, and data 
recorder/instrumentation operator, plus all the equipment, 
was running out of room! Also, trying to monitor and record 
simultaneously, with minimal acceptable accuracy, multiple- 
point thermcouple outputs on two different instruments was 
difficult. 

The test was viewed as a success despite the problems 
with the core box, towing speed, and crack development. A 
larger electric motor was used in future experiments to 
increase the power available to overcome excess drag. The 
freezing technique reverted to the original 1-2 inch layer 
thicknesses. If additional wake temperature data were to 
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be obtained, the pontoon bridge would require modification 
so that it would remain exactly behind the ice model and 
not drift to either side. The pontoon vehicle did not 
follow in the wake of the ice model. In addition, a second 
portable strip chart recorder was obtained to reduce the 
switching complexity when recording multiple-point boundary 
layer temperature profiles. 

6. Test No. 6, April 27, 1978 

In addition to the overall test objectives of previous 
experiments, two additional areas of investigation were 
considered during this test. The first area was the 
determination of the effect, if any, of the turbulence 
created by the towing cable and thermcouple wire bundle on 
the regression rate of the bow of the model. The second 
area was the examination of the difference in melt rates at 
two points on a side ripple. 

The ice was frozen using the 1-2 inch layer thickness 
technique with the core box and thermcouples placed as 
shown in Figure 60. The box had been repaired from the 
damage incurred in Test No. 5, and the towing bollard was 
added as an integral part of the structure of the box. The 
model appeared in excellent condition at the launch site as 
illustrated in Figure 61, which also shows the thermcouple 
positions for this test. The channels were assigned as 
follows: stagnation point boundary layer (11), bottom 

surface boundry layer (12) , ice center forward (13) , forward 
ripple boundary layer (14 and 15) , rear ripple boundary 
layer (16 and 17) , and rear ice center (18) . 



Previous ripple data measurements indicated the 
average peak-to-peak wavelength of the ripples was approxi- 
mately 6 inches. The melting rate varies at different 
points between the crest and trough of a ripple. The arrange- 
ment of thermcouples was an attempt to measure the different 
melt rates due to a ripple. Two stations were selected 
where boundary layer profiles previously had been recorded. 

A pair of thermcouples was positioned at each station with 
a 4-inch separation. Initially the sensory tips were at 
identical distances from the receding edge and at equal 
depths in the ice. Both thermcouples were monitored simul- 
taneously on a real-time basis by a two-channel strip chart 
recorder. Since the separation was less than one ripple 
wave length, the position of each tip relative to the ripple 
was indicated by the sequence of exposure and shape of the 
temperature profile. 

The thermocouples were placed in the ice to achieve 
the following exposure sequence: stagnation point, forward 

side station, rear side station, and bottom station. Test 
No. 6 was the most ambitious attempt at multipoint boundary 
layer recording throughout the series of tests. However, 
launching difficulties during the test diminished the thermo- 
couple data. The difficulties began when the ice container 
was placed in the water and prematurely lowered before the 
ice could free itself. This resulted in two adverse conse- 
quences; first, the thermocouple electrical connections were 
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flooded effectively nullifying their output, and, second, 
the ice separated from the container unevenly stern-first 
putting excessive stress on the forward section; as a 
result, the forward section of ice broke at a location 
behind the forward ice center thermocouple as pictured in 
Figures 62 and 63. The decision was made immediately to 

disconnect the thermcouple wires from the forward section 

and continue the test as planned with the remaining section 
of the model. 

The ice model was pushed from the rear rather than 
towed during this test, as pictured in Figures 64-66, to 
determine the effect of the towing cable and thermcouple 
wire bundle on the development of the bow erosion pattern. 
The wood pushing bar was hinged on either end to allow 
vertical movement of the model and boat. The bar was bolted 
with a facing plate to the rear bollard while clamped to 
the bow of the boat. Two torque lines, running from both 

arms of the main bollard back to the line handler on the 

boat, aided in directional control of the model and mini- 
mized the torque applied on the rear bollard and bow of the 
boat. The underway portion of the test was continued until 
significant undercutting was observed on the model, at which 
time the ice was recovered as shown in Figures 67 and 68. 

Observations of the model on the pier showed that 
the same bow surface erosion characteristics developed to 
the same extent on this test as when the model was towed. 
Although the bow erosion did not have exactly the identical 
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shape because of its slanted front as a result of the break, 
close examination showed clearly that the same process of 
erosion had taken place. The well-defined periodic ripples 
were present as in previous tests. An improved method of 
measuring these ripples was utilized as shown in Figures 69- 
71. A wooden reference frame was prefabricated and posi- 
tioned in slots on the core box. A Longfellow Linear Motion 
Position Transducer Model LF-6-101, pictured in Figure 70, 
was mounted on either rail of the frame and moved the entire 
length of the frame is illustrated in Figure 71. As the 
transducer moved along the rail, the position arm moved 
over the ice surface, and the resulting changes in trans- 
ducer output voltage were recorded by a strip chart recorder. 
The end of the chart paper was attached to the transducer 
chassis so that the transducer moved along the rail at the 
same speed as the chart paper moved through the recorder. 

Thus the recorder printed an exact geometrical profile of 
the ripples over the length of the model. 

This test was the last in the series to study a 
single ice model. The knowledge and experience gained from 
the series were applied to another test in which two similar 
models were towed in tandem to determine the effect of the 
first model on the melting rate and surface regression 
characteristics of the second. This test and the work 
leading to it are discussed in Chapter V. 
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F . COMMENTS 



An invaluable store of knowledge and experience on 
testing large scale ice models underway in seawater was 
gained as a result of the series of experiments in Monterey 
Bay. While not every test objective was met on every test, 
those that were accomplished added to the experience that 
was nonexistent prior to this series. 

One important point was highlighted by the test series. 
That was the necessity of controlled laboratory conditions 
for experiments with the wide scope of the tests. The 
arrangements alone for the ice model, transportation, launch 
site, boats, equipment, and assistance required the coordina- 
tion of six different agencies. Add to this the unpredic- 
tability of weather, and one soon grasps the extensiveness 
in preparing for a single test. The complexity involved 
during the tests in monitoring instrumentation while con- 
trolling the boat and model has been previously discussed. 
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V. TANDEM ICE MODEL EXPERIMENTS 



A . BACKGROUND 

The ice model test series described in Chapter IV showed 
that the most intensive melting and subsequent regression 
occurs at the bow of the model. This rapid melting is a 
function of pressure gradients due to variation in the local 
velocity and the higher local temperature gradient profile 
on the bow of the model. The "cavern"-shaped eroded area 
on the lower frontal surface of the model is that area at 
which the higher local velocity and acute pressure gradients 
are most pronounced, which results in the high melting and 
regression rate in this portion of the bow. Griffin's 
[1,2] model showed the strong dependence of melting rate 
on the differential between the ice surface temperature and 
the local ambient seawater temperature. 

One possible solition to minimizing the pressure and 
velocity gradients and the ice surface-seawater temperature 
differential is to tow the iceberg in the wake of a 
"sacrif icial'' iceberg. The sacrificial iceberg could shield 
the rear iceberg by lowering the velocity of the seawater 
impinging on the bow and also provide a plume of seawater 
cooled to a lower temperature than the ambient seawater 
temperature. The former would reduce the pressure gradi- 
ents on the bow while the latter would reduce the ice surface- 
seawater temperature differential, possibly resulting in a 
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substantial reduction of the high regression rate of the 
rear iceberg. 

While the general theoretical heat transfer principles 
were well-known and understood, previous research in this 
area with ice bodies underway in seawater was nonexistent. 

As a concluding experiment for this thesis and an introduc- 
tion into possible future work, a test was conducted in 
Monterey Bay consisting of towing two large-scale ice 
models in tandem on December 10, 1978. The models were an 
identical pair and similar in size and construction to those 
in the previous single-model series. 

B. LABORATORY TESTS 

Due to limitations on time and resources, only one large- 
scale experiment on tandem-towed models could be conducted 
for this thesis. In order to obtain decisive results from 
the large-scale test, a series of laboratory experiments 
were conducted first on small-scale models utilizing the 
water-channel pictured in Figure 5. The characteristics 
of this water channel are listed in Table I. The primary 
issue to be addressed in these experiments was the distance 
the second model should be towed behind the first to illus- 
trate most drmatically the "protection" the bow of the 
second would derive from the presence of the first. Of 
practical concern was the realization that full-size ice 
bodies could not be towed so close in tandem that longi- 
tudinal variations in their individual velocities could cause 
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an "accordian" effect which would result in catastrophic 
collisions between the icebergs. 

Froude number scaling was utilized to model the Monterey- 
Bay ice model for the laboratory experiments. In Chapter III 
it was stated that for Froude number scaling: 



Fn 



model 



Fn 



prototype 



while eq. (2) showed that 



( 9 ) 



/gL 

with velocity (V) and model length (L) as the controlling 
parameters. Substituting: 



(-=) 

/gL model 



/gL prototype 



( 10 ) 



rearranging ; 



V model ^prototype ^L, 



mode 1 



1/2 



( 11 ) 



prototype 



The length and velocity of the prototype (16 ft., 1.78 ft/sec) 
are known. A freezing container was available with dimen- 
sions of L = 13.5 inches, W = 4.5 inches, and D = 2.5 inches 
for forming the laboratory model. From eq. (11) the model 
velocity was 0.472 ft/sec. The model Froude number from 
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eq. (2) was 0.023, thus precluding attached bow waves. The 
low channel velocity eliminated any possibility of signi- 
ficant blockage or wall interference. Although the width 
and depth of the model was not scaled exactly relative to 
the length, their sizes were satisfactory for the laboratory 
tests since drag of the model was not at issue. 

The laboratory models were frozen in a standard freezing 
compartment for a two day period at a temperature of 0°F. 

The model was frozen as a single 2.5-inch layer with two 
wooden handling devices. Two devices were frozen in each 
model to assure lateral directional stability during the 
tests . 

The laboratory tests were conducted by placing the ice 
models in the flow channel and allowing them to float using 
their own natural buoyancy as pictured in Figure 72. Direc- 
tional control was attained by forcing the wooden handling 
devices to ride in a slot formed by two guide rails that ran 
longitudinally the length of the channel centerline. Each 
test was conducted for a continuous period of two minutes. 
This two-minute period was determined by placing a single 
model in the channel for a sufficient period to allow 
decisive erosion of the bow as pictured in Figure 73. The 
rear model was placed initially at a separation distance 
corresponding to one model width (W = 1.5 inches), the 
characteristic model dimension chosen for the tests. On 
each successive test the rear model was moved incrementally 
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closer reducing the separation distance. During each test 
the flow immediately adjacent to the front and rear models 
was studied by injecting a colored dye in the water in close 
proximity to the area of interest as shown in Figure 73. 

C. RESULTS OF LABORATORY TESTS 

The decisive "necking" observed in the initial laboratory 
experiments was also noted during this series of experiments. 
This "necking" effect served a very useful purpose during 
the tandem laboratory experiments by providing a well-defined 
bow erosion signature for comparing erosion on the two models. 
Chapter I and Chapter VI describe the "necking" on the 
laboratory models and the flow characteristics which cause 
this rapid erosion to occur. The flow recirculation was 
seen easily with the aid of the injected dye as pictured 
in Figure 73. 

When placed at the initial separation distance of one 
model-width (W) , the rear model bow eroded to exactly the 
same shape as the forward model bow. The rear model appeared 
to have derived no substantial benefit from the presence of 
the forward model. At separation distances of 2W/3 and W/2 
the rear model bow eroded in the identical shape of the 
forward model, but at a slower rate. At a separation dis- 
tance of W/3, the rear model bow eroded in a decisively 
different manner from that of the forward model as shown in 
the sequence of Figures 74 -76. The decisiveness of this 
change in erosion pattern can be observed further by examining 
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Figures 77 and 78, which show the bottom surfaces of the 
forward and rear model respectively. The complex flow 
pattern between the two models at the W/3 separation distance 
is shown in Figure 79. 

As a result of this experiment in the laboratory, the 
large-scale tandem-towed model test in Monterey Bay was 
conducted utilizing the W/3 separation distance between 
models. Shorter separation distances were not examined due 
to the practical uncertainties of longitudinal oscillations 
with resulting collisions mentioned earlier in this chapter. 

D. MONTEREY BAY TANDEM TEST, DECEMBER 10, 1978 

The models used in the tandem experiment were frozen 
in identical fashion, utilizing the 1-2 inch layer freezing 
technique with similar core boxes installed. Thermocouples 
were not installed in the ice due to the complexity of the 
experiments in handling two large ice models simultaneously. 
Since only one crane was available, the models were launched 
separately (forward model first) and then linked together 
in the Bay as pictured in Figure 80. The linking was accom- 
plished by two rigid metal rods connected to either arm 
of the towing bollards that were sized to provide a constant 
W/3 (16 inch) separation distance between the models. 

The models were towed over the same course as in previous 
experiments utilizing the same Boston Whaler towing vessel 
as pictured in Figure 81. While the rear model remained 
behind the lead model for some brief periods (Figure 82) , 
lateral instabilities in the models resulted in lateral 
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motions and misalignment of the models during the major 
portion of the experiment as shown in Figures 83 and 84. 
Attempts at restoring the misalignment by providing a drag 
force from a following Boston Whaler on the model train 
(Figure 84) were unsuccessful. 

Recovery was initiated after 21 minutes underway when 
significant undercutting was observed on the lead ice model. 
The rear model was recovered first, followed by the lead 
model. Both models were again weighed to determine total 
net weight lost during the test. 

E. RESULTS OF TANDEM BAY TEST 

While decisive indications were obtained from the 
experiment, the results were in a large measure inconclu- 
sive due to the nature of the experiment. Weight indications 
showed that the lead model lost 2750 lbs. and the rear model 
2000 lbs. While this large weight loss differential would 
appear to prove decisively the "sacrificial" iceberg theory, 
one must remember that the lead model was in the water an 
additional 36 minutes due to the limited single launch and 
recovery capability available. 

The recovered lead and rear models are shown in their 
bow view in Figures 85 and 86 respectively. The bow erosion 
of the rear model is approximately the same as the lead 
model which would negate the shielding theory provided by 
the lead; however, the rear model was underway eight addi- 
tional minutes over the lead in order to link up with the 
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naturally drifting lead model in the Bay. Additionally, 
the misalignment resulted in the rear model's bow being 
exposed to the full flow of the oncoming seawater during the 
major portion of the test. 

The Tandem Bay test procedure can be improved by 
utilizing two lifting cranes to effect a simultaneous 
launch of both models and by connecting the tandem models 
through a lightweight rigid frame to prevent lateral motion 
between the models. The former would be more expensive and 
require additional pier space, but the result would be a 
more accurate reflection of the actual underway mass loss 
differential between the models. The latter would assist 
in better defining the bow erosion differential between the 
models . 
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VI. EXPERIMENTAL RESULTS WITH THEORETICAL CORRELATIONS 



A. GENERAL 

The following are discussions of the experimental 
results determined from the ice model test series described 
in Chapters III-VI. Comparisons with appropriate theoretical 
predictions or other experimental results are made including 
descriptions of the theoretical models and experimental 
techniques . 

B. COMPARISON OF EXPERIMENTALLY MEASURED ABLATION RATES 

WITH THEORETICAL PREDICTIONS 

Griffin developed a model for predicting the ablation 
rates of icebergs moving through seawater in laminar flow 
[1] . He modified his model [2] to account for turbulent 
flow in conjunction with the large-scale ice model experiments 
conducted at NPS. 

The analysis model considers turbulent flow over a 
flat plate with zero pressure gradient and solves the boun- 
dary layer equations for simultaneous heat, mass, and 
momentum transfer. The two-dimensional problem is reduced 
to one dimension by integrating across the boundary layer. 

The phase transformation from solid to liquid occurs under 
steady-state conditions enabling a coordinate system to be 
affixed to the melt interface. A constant density is 
assumed, and the ice is a pure solid phase. Gravity is 
neglected, thus neglecting the buoyancy effects on the heat 
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transfer. The model requires inputs for the physical 
properties of seawater and ice, including the applicable 
dimensions and free-stream velocity. The properties of 
ice and seawater are assumed to be constant for the solution. 

The graphical results for overall surface regression 
for three of the large-scale experiments in Monterey Bay 
are shown in Figures 87 - 89 and summarized in Table V. 

A comparison of the actual regression rates at three points 
on the model tested on October 2, 1977, with the prediction 
from Griffin's Computer Code is shown in Figure 90. The 
agreement between experimental values and theoretical pre- 
dictions is good. A plot of turbulent flow Stanton numbers 
versus Reynolds numbers is shown in Figure 91 for the 
range of Prandtl numbers typical for seawater. Data points 
from Table V of NPS experiments are shown again to be in 
good agreement with the theoretical curves. A similar 
plot of Stanton numbers versus Prandtl numbers is shown in 
Figure 92 for a range of Reynolds numbers from laminar to 
turbulent flow conditions . 

The conclusion drawn from the above comparisons was 
that the Griffin turbulent-flow iceberg ablation model was 
satisfactory for predicting the ablation rate at various 
points on the sides an iceberg given the appropriate physi- 
cal conditions such as ice and water properties, velocity, 
and water temperature. However, the bow area of the iceberg 
requires further treatment. 
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STAN.T0N NUMBER = HEAT TRANSFER 00EFFICIENT/(IENSITY X SPECIFIC HEAT x TOW SPEED) 
* Values for Seawater 



C. MEASUREMENT AND VERIFICATION OF THERMAL BOUNDARY 

LAYER PROFILE 

The profiles of the water temperature versus distance 
from the ablating surface for experiment Nos. 4 and 5 in 
Monterey Bay are shown in Figures 93-95. The data on 
Experiment No. 4 were obtained by recording a digital readout 
from the thermocouple array described in Chapter IV while 
the data for Experiment No. 5 were recorded by a single 
thermocouple output to a continuous strip chart recorder. 

Due to the integral solution technique applied across 
the mass diffusion, thermal, and momentum boundary layers, 
Griffin's model does not provide predictions for the 
thermal boundary layer profile. However, an approximate 
determination of the accuracy of the initial profile slope 
could be calculated by applying a one-dimensional energy 
balance across the ice-water interface as shown in Figure 



107. 



From Figure 107: 



heat transfer from water ( 




heat transfer into solid ice 



cal 
Sec cm 



M 



O 3. X 

latent heat of fusion ( 



The energy balance is : 



p . V M + q . 
is ^1 



( 12 ) 
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where 



ice density 




V g = regression rate (^-) 



Also : 



■‘w 



= K (il) 

w ( 3y j w 



(13) 



and 



= K. 

l 



(£ 1 ) . 

v 3y ; i 



(14) 



where K is the respective thermal conductivity and ^ — is the 

dy 

respective initial slope of the temperature profile. 

From Ref. [7] at 0°C: 



K. = 5. 35 x 10 

l 



-3 



cal 

’C cm-sec 



K = 5.33 x 10 

W 



-3 



cal 

’C Cm-sec 



M = 79.77 



cal 



p = .9168 g/cm' 
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The ice interior temperature Profile - Test No. 4: 



From Figure 96; 



(£1) . 
l 3y ; i 



1.1 °c 



(.0625 in) (2.54 

in 



From Figure 94; 



(IT) 
v dy vi 



6 °C 



(.04 in) (2.54 ™) 
m 



From Table V; 



V s = 5.9 in/hr (at L = 2.4 rn) 



Solving eq. (14) and substituting into eq. (12) gives 



A , . cal 
%, ‘ °' 34 1 



sec-cm 



Solving eq. (13) gives a value of 



q 



w 



0.3 



cal 

7 ' 

sec-cm 



Reasonable agreement occurs between theoretical and observed 
results indicating that the measured profile is a reasonably 
accurate representation in at least the critical initial 
slope . 
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D. RIPPLE EFFECTS AND MEASUREMENT 



1 . General 

Throughout the test series, ripples were observed 
to have formed on the sides of the ice models tested in 
Monterey Bay. The ripple formation is not unusual as turbu- 
lent flows in nature are observed to produce several differ- 
ent types of regular waves, including wind-generated water 
waves. Carey [8,9], Larsen [10], Ashton and Kennedy [11,12], 
and Hsu [13] have described and analyzed the occurrence and 
behavior of regular waves on the underside of river ice 
covers. These waves are known as ice ripples and can be 
expected to form at any ice surface exposed to a turbulent 
water flow, at least when the boundary is receding due to 
melting. Tatinclaux and Kennedy [14] developed an analyti- 
cal model for formation of ripples on icebergs. 

2 . Ripple Effects on Drag 

Both heat and momentum are transferred between the 
ice-water boundary layer and the surrounding flow by means 
of turbulent exchange processes. In addition, the ripples 
eventually will become sufficiently high so that the flow 
will separate near the crest of each ripple. The normal 
relatively low pressure at the crest will then continue 
across the rear face of the ripple until the flow reattaches 
at the frontal face of the next ripple. The large pressure 
differential across the ripple results in a significant 
drag force exerted on each ripple. Also there will be 
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gravity waves generated in the water by the ripples, thus 
adding to the iceberg drag. 

Several methods have been proposed for calculating 
the composite friction coefficient by Carey [8] , Larsen 
[10] , and Uzuner , et. al. [15] . Larsen showed that drag 
increases almost linearly with ripple steepness (peak-to- 
peak ripple height divided by ripple length) up to a value 
of about 0.1 and much more rapidly thereafter. For fully 
developed ice ripples, values of steepness in a range of 
0. -0.15 are not uncommon. Further, such ripples could 
increase frictional resistance by 50 - 100 percent. 

3 . Ripple Effect on Heat Transfer 

The formation of ripples at the ice boundary has 
the effect of increasing the rate of heat transfer at the 
boundary. The increase in heat transfer is due, first, to 
the increased hydraulic roughness which affords a larger 
net area available for heat transfer and, second, to the 
increased turbulent intensity and flow separation caused by 
the ripples. 

4 . Measurement of Ripples on NPS Tests 

Chapter IV described the various tehcniques used to 
measure the ripples on the sides of the ice models tested 
in Monterey Bay. From the data gathered, histograms of 
ice ripple wavelength were constructed and are shown in 
Figures 97-99. From the histograms the average wavelength 
was determined for each experiment. 
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Experimental programs conducted at the Iowa Insti- 
tute of Hydraulic Research on the formation of ice ripples 
have been reported by Ashton and Kennedy [12] and Hsu [13] . 
The experiments were conducted on an ice slab of dimensions 
L = 40 ft. and W = 2 ft. over a velocity range of 0.30 ft/sec 
to 2.7 ft/sec. Hsu reported that the average wavelength 
(X) is inversely proportional to the mean flow velocity 
(U) and that the Reynolds number based on wavelength is 
constant and given by: 

n \ d 

R = — = 6.2 x 10 (15) 

v 

The standard deviation of the experimental results from this 

4 

average is 1.1 x 10 . 

The average wavelength for each NPS experiment was 
determined from eq. (15) and compared to the corresponding 
average wavelength determined from the histograms. The 
results are listed in Table VI and show good agreement. 

E. THEORETICAL PREDICTIONS AND MEASUREMENT OF ICE 
INTERIOR TEMPERATURE PROFILES 

The temperature profile within the ice model was 
measured during the underway portion of Test No. 4 in 
Monterey Bay. The thermocouple array discussed in Chapter IV 
and shown in Figures 41 and 42, was read at specific time 
intervals for a period in excess of 37 minutes; readings 
were voice recorded. 
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TABLE VI 



COMPARISON OF RIPPLE DATA MEASUREMENT 



DATE: 



10/3/77 12/10/77 



Velocity (U) 1.5 

(ft/sec) 



1.4 



Wavelength 5 . 4 

U 1 ) (in) 



6.1 



Wavelength 5.4 

U 2 ) (in) 



5.8 



Note : 

Wavelength from eq. (15) 

A 2 Wavelength from histograms. 



4/27/78 

1.9 

4.4 

4.7 
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The transient thermal conduction equation for one 
spatial dimension was solved, and the results were used to 
correlate the experimental measurements . The equation and 
development of the solution is given in Appendix B. Included 
in Appendix B is a Hewlett-Packard-9831-A computer program 
to generate the theoretical temperature profiles, thermal 
diffusivity of ice (a) , and ice regression rate (V g ) . 

The comparison of the theoretical and experimental 
temperature profiles is shown in Figures 100 - 105 for each 
respective time interval. The theoretical curves represent 
a best parametric fit with the experimental values by 
independently varying such parameters as thermal diffu- 
sivity and regression rate as described in Appendix B. The 

resulting theoretical thermal diffusivity and ice regression 

-7 2 

rate from the best fit condition were 3.0 x 10 m /sec and 

-5 

2.2 x 10 m/sec respectively. The actual thermal diffu- 

-7 2 

sivity of ice determined from Ref. [7] was 5.3 x 10 m /sec 
for pure solid ice. The actual regression rate, determined 

by initial and final ice model dimensions for the 37 minutes 

-5 

of underway time, was 3.53 x 10 m/sec. 

Figures 100 - 105 show compatability between the theoreti- 
cal and experimental data for temperature profiles within the 
ice. The theoretical curves and experimental values do 
diverge to some extent at time intervals 5 and 6. One 
hypothesis for the divergence is the significant change in 
thermal diffusivity in the relatively steep temperature 
gradient. 
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The variance between theoretical and actual ice thermal 
diffusivity is explained by Mellor [16] as the result of 
air bubble entrapped in the ice which lowers the thermal 
diffusivity from that of pure bubble-free' ice. There was 
a significant bubble population throughout each model tested 
in the bay. 

Since the time over which the regression rate was calcu- 
lated did not include the 7 minutes the ice was motionless 
in the water, it was not surprising that the measured 
regression rate was higher than the theoretical prediction. 

F. MEASUREMENT OF ICE MODEL WAKE TEMPERATURE 

The method used to measure the wake plume on Test No. 5 
in Monterey Bay is described in Chapter IV. The method 
proved largely unsuccessful for the following reasons. 

First, the pontoon bridge did not remain aligned continu- 
ously with the centerline of the model throughout the test 
which limited the opportunities to obtain data. Second, the 
thermocouples had to be read individually in stepwise fashion 
due to recording equipment limitations thus adding variance 
to the recorded temperature profile at any given time. 

Third, variations in surface water temperatures, large 
scale eddies in the wake, wave motion and similar effects 
added to scatter in the measured values. The recorded 
temperature profiles are shown in Figure 106 for five 
different time intervals. A search of the literature showed 
no data for similar measurements for comparison with the 
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NPS experimental results. Further, efforts at measuring 
the wake plume were not attempted due to limited time and 
resources . 

The conclusion is that measurement of wake temperatures 
is not as straightforward as one might expect. While the 
equipment recording problem could be solved by continuous 
monitoring and recording techniques, the large scale eddies, 
which were observed in the wake with dye injection, signi- 
ficantly complicate the interpretation of temperature 
profiles obtained. 

One obvious reason for measuring the ice model wake 
temperatures is to aid in the determination of the effective 
cold-water "insulation" provided by the model for a similar 
model being towed in the wake. 

One alternative method for measuring the wake tempera- 
ture would be to measure it with the second model in tow. 

The close proximity of the bow of the trailing model to 
the stern of the lead model would have a stabilizing effect 
on the flow leaving the lead model. An attempt at this 
method was planned for the tandem test in Monterey Bay but 
could not be executed because of the lateral oscillations 
of the models. 

G. FLOW SEPARATION ON THE ICE MODELS 

The decisive "necking" observed on the small-scale 
laboratory ice models was due to the recirculating flow in 
the forward one-third region on the models. The recirculation 
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was a result of the initial flow separation around the 
sharp frontal corners and subsequent reattachment as des- 
cribed in Chapter III. The separation region, where the 
necking occurred, was observed visually with use of dyes 
as described in Chapter V. However, the "necking" was never 
observed on the large-scale ice models tested in Monterey 
Bay, despite the fact that flow separation around the initial 
90-degree corners certainly occurred. A discussion of the 
flow separation conditions for the two modeling situations 
follows . 

Consider first the frontal corners on the small- and 
large-scale models. The corners are initially areas of 
intense heat transfer and erode rapidly once the model starts 
moving in the fluid medium. As the corners erode, their 
radius of curvature increases thus decreasing the tendency 
for flow separation. The actual rate of decrease in the 
flow separation tendency will vary depending on the laminar 
or turbulent condition of the flow. In the laboratory tests, 
the flow was laminar due to the low Reynolds numbers involved. 
The small frontal area of the laboratory model did not offer 
sufficient opportunity for the flow to reach turbulent con- 
ditions prior to rounding the corners. Laminar flow is 
more sensitive to flow separation even on relatively stream- 
lined bodies. Therefore, although the corners became 
rounded on the small-scale model, the flow, which remained 
laminar, continued to separate with the subsequent "necking" 
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as previously seen in Figures 6 and 7 of Chapter III and 
Figure 73 of Chapter V. 

On the large-scale model, the flow was turbulent as 
evidenced by the full-length ripple system on the sides and 
"cavern"-shaped area on the bow of the model. Due to the 
significantly larger frontal area on the large-scale model 
as compared to the small-scale model, the flow has an 
enhanced opportunity to reach a turbulent condition prior 
to rounding the corners. A well-known fact is that turbu- 
lent flow has a lower tendency than laminar flow for separa- 
tion. Therefore, one hypothesis is that once the corners 
on the large-scale model become rounded, the flow separation 
ceased and no observable "necking" occurred. Such a hypothesis 
could be tested using both potential flow theory and actual 
continuous observation of the large-scale model as was done 
in the laboratory experiments. 

The significant difference in the melting shapes of the 
small- and large-scale models emphasizes the importance of 
careful treatment of the scaling factors such as Reynolds 
number when modeling a large iceberg. The melting shapes 

were obviously strongly dependent on Reynolds number, which 

4 6 

was on the order of 10 in the laboratory tests and 10 

in the Bay tests. 
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VII. RECOMMENDATIONS FOR FUTURE WORK 



A. GENERAL 

There are many areas of research concerning glacial 
ice moving through seawater which strongly merit future 
examination. Included are improvements in testing techniques 
and instrumentation used to obtain such basic data as drag 
force, velocity, temperatures, and salinity. The following 
are brief discussions concerning these research areas. 

B. TOPICS FOR FUTURE RESEARCH 
1 . Measurement of Drag 

For a square bow, such as shown in Figures 27, 29, 
and 31, the drag is dominated by form drag. At the low 
Froude numbers, the wave drag is negligible. Skin friction 
is relatively insignificant. 

As the iceberg melts, the shape tends to a geometry 
yielding lower form drag. Also the bow may be shaped to 
give lower form drag. In that case, skin friction becomes 
more significant. There is an obvious interaction between 
convective heat transfer and skin friction. One expression 
of the relation is Reynold's analogy. 

As a result of the past towing experience, the NPS 
team should develop techniques to measure drag accurately. 
Drag information is useful in the interpretaion of shape 
changes and for correlation with melting rates. 
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2 . Tandem Towing 



Only one test was conducted at NPS in this inter- 
esting and complex method of moving icebergs. One potential 
advantage of tandem towing is the fluid mechanical insulation 
(cold water) provided by the lead iceberg. The lead iceberg 
can be considered to be a "sacrificial" iceberg providing 
the cold water for insulation. Evaluation should include 
boundary layer temperature profile measurement and flow 
field surveys at the rear of the lead iceberg. 

3 . Influence of Buoyancy on Heat Transfer 

The paper by Josperger [17] clearly demonstrates 
the importance of buoyancy forces on natural convective 
heat transfer. The critical Grashof number for transition to 
turbulent flow is a factor of 10 lower for fresh water ice 
melting in seawater than for a vertical hot plate in air. 

An analytical model should be developed to determine 
the influence of buoyancy. A buoyancy term should be 
introduced in the equations of motion developed by Griffin. 
The impact of this term on heat transfer should be examined. 

The Naval Postgraduate School is about 60 miles from 
San Antonio Reservoir which is near King City, California. 

In order to test the theory concerning buoyancy and heat 
transfer under forced convection (i.e., moving iceberg), 
several icebergs could be towed in San Antonio Reservoir. 

A comparison could be made between melting rates in fresh 
water and seawater. Furthermore, the inverse melting 
problem, i.e., an iceberg of frozen seawater melting in 
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fresh water, can be studied to gain insight on buoyancy 
effects in forced convection. 

4 . Instabilities of Insulating Material, e.g., 

Flagwaving Instability 

Consider an insulation scheme using thin plastic 
sheets. A protective layer of cold water develops between 
the iceberg and the plastic sheet. 

There are several excellent books , e.g., Lin [181, 
Chandrasekhar [19], and Shen [20], which discuss stability 
of flows. In rotating coordinates, the book by Greenspan 
[21] is especially relevant. One well-known instability 
is that which occurs in parallel flows of fluids with 
unequal densities. 

A survey is needed of the instabilities that may 
arise in the context of icebergs with insulation. Such 
a survey is proposed here. 

5 . Solid Mechanics of Icebergs 

The Naval Postgraduate School has a Chair in Arctic 
Research. The incumbent starting in July, 1978, is Dr. 
Wilford F. Weeks, who is a Research Glaciologist with the 
U.S. Army Cold Region Research and Engineering Laboratory. 
Dr. Weeks has agreed to devote some of his time to the 
question of solid mechanics of icebergs. 

It is apparent from the technical papers at the 
Iceberg Utilization Conference that many of the existing 
engineering techniques relative to cracks, fractures, stress 
distribution, multiaxial stresses, etc., 
iceberg structures . 



are applicable to 



6 . Development of Instrumentation 



To verify the computer code of Griffin, measurement 
of temperature, salinity, and velocity profiles would be 
especially helpful. 

Velocity of the iceberg is measured using standard 
oceanographic equipment. Drag is measured also using 
standard equipment. 

Measurement of salinity and velocity profiles 
requires special consideration. The surface of the ice is 
receding. The iceberg is a remote vehicle. 

As a research task, techniques should be developed, 
or improved, for measuring: 

1. velocity of iceberg 

2 . drag 

3. salinity profile 

4. velocity profile 

5. temperature profile 

As a further aid to the experiment, underwater photographs 
and underwater movies of an iceberg in motion should be 
made . 

7 . Laboratory Tests 

Towing tests in Monterey Bay are time consuming and 
expensive. Hence, each test must be planned carefully and 
overall test objectives must be defined clearly. To supple- 
ment the towing tests in Monterey Bay, tests can be conducted 
using the water tunnel and water channel located in the 
Department of Mechanical Engineering. 



By using a larger test facility, larger test 



specimens, and higher flow velocity, turbulent flow may 
be achieved in the laboratory. Special emphasis should be 
placed on defining the usefulness of boundary layer trips. 

Due to the relative simplicity of laboratory testing, 
many tests can be conducted in the laboratory compared to 
larger scale towing tests in the ocean. 

8 . Field Tests 

There have been three options for large scale tests 
discussed and evaluated by the Naval Postgraduate School 
research team. Option 1, which is towing in Monterey Bay, 
has been exercised six times. Option 2 involves use of 
tidal currents at Moss Landing. Elkhorn Slough has a 
narrow opening; the bridge for State Highway 1 crosses the 
Slough at this point. Due to tides, a velocity of 1.0 knot 
occurs and persists for about 30 minutes. Using the bridge 
as a point to tie the two ropes, tests could be conducted 
with the 16-foot iceberg. The bridge and the nearby 
navigational pilings could be used as observation points. 

Option 3 is to tow the large scale iceberg in San 
Antonio Reservoir. The goal of these tests would be to 
determine influence of buoyancy on heat transfer. Also at 
San Antonio Reservoir water temperature would allow use of 
swimmers near the iceberg; the swimmers could observe the 
flow at close hand. 

Experience has shown that the field tests would be 
facilitated greatly if an observer could ride the iceberg. 
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Hence a double wall box is proposed. For the size of manned 
iceberg, the weight of the ice is 7735 lbs. With a payload 
of 1500 lbs, which includes the weight of the inner box, 
the freeboard is almost 12 inches. The design uses a broad 
beam to avoid stability problems if ice were to break away. 

Advantages of a manned iceberg include the following: 

1. ability to observe the flow. 

2. ability to conduct experiments, start recorders, 
and vary profile instruments. 

9 . Insulation 

Many different forms of insulation have been proposed 
including plastic sheets, frozen wood chips, plastic foam, 
etc. An alternate scheme might be the use of air bubbles 
at locations with high heat transfer rates. 

To examine the concept, an experimental program is 
proposed. First, the bubble insulation would be investigated 
in the water channel. Second, one or more of the large 
scale icebergs would be dedicated to an evaluation of 
bubbles for insulation. Third, depending on the outcome 
of the experiments , a study would be made of pumping require- 
ments. Also, the study would examine installation of 
bubble source on the iceberg. 

10 . Potential Flow Computer Code 

Pressure gradients have a pronounced effect on boun- 
dary layers. A potential flow computer code would be used 
to calculate pressure coefficients and pressure gradients. 
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9 



The shape of the body will be that shape observed during 
experiments conducted in Monterey Bay. 

A computer code has been developed by Professor 
C.J. Garrison under a NSF grant. Professor Garrison is a 
faculty member in the Mechanical Engineering Department at 
NPS . 
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APPENDIX A 



AREAS OF CONCERN WITH ICEBERG DETERIORATION 

The following is a discussion by Stolfi, et al [22] of 
the motivation and scope of interest in iceberg deterioration. 

A large percentage of the land surface of the globe 
consists of arid and semiarid regions which support less 
agriculture and smaller populations than other areas. In 
certain semiarid regions where agriculture is successful, 
that activity is dependent upon well-developed irrigation 
technology and the light variable rainfall of the region. 

The arid regions, and especially the areas in them lying 
far from the world's oceans, present extraordinarily challenging 
problems for development in some future time frame. Certain 
arid regions and many semiarid regions which lie close to 
the world's oceans contain significant populations and 
present special- opportunities for expanded development and 
insurance against the periodical, economically damaging 
droughts associated with rainfall scarce areas in which 
productive agricultural activity is being conducted. 

Several appraoches can be envisaged for the support of 
agriculture and associated populations and industry in arid 
and semiarid regions. The light, existing rainfall can be 
gathered, stored, managed, and metered out to the surrounding 
countryside and cities inwater storage and irrigation works 
which conserve and distribute the existing rainfall. Similar 
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water storage and distribution works can transport water 
from heavier rainfall regions located economically close 
of the drier areas, to make the water available for agricul- 
ture and maintain the water tables and other underground 
water resources. In cases of deep and extensive drought, 
however, as in Western North America, 1975-1977, the lack of 
rainfall will empty reservoirs, lower permanent rivers, and 
drive down water tables over large regions threatening 
immense damage to grazing, agriculture, and ultimately 
industry. In both arid and semiarid regions in which 
productive agriculture has been developed through water 
conservation, storage, and distribution technology, the 
water supplies can be supplemented by (1) inter-region trans- 
fer of water from distanct heavier rainfall areas, (2) 
fossil-fuel, nuclear, or solar desalinization of ocean 
water, and (3) transfer and terminal processing of Antarctic 
ice. The approaches can be summarized as follows: 

Support of Agriculture and Industry 
in Water Deficit Regions 



Arid Regions Semiarid Regions 



(Sahara, Namib, Atacana, etc.) (California, Argentina, South Africa, etc.) 



Transfer of Water frcm Closely Conservation, Storage, Distribution 

Adjacent Areas with Existing Rainfall of Light Rainfall Resources 



Inter-Region Transfer Ocean Desalinization Antarctic Ice 

of Whter frcm Distant (Using Fossil, Nuclear Ocean Transport 

Areas or Solar Energy) and Terminal Processing 
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The expansion of economic activity into arid regions and 
the intensification and assurance of uninterrupted activity 
in semiarid regions demand advances in technology related 
with Inter-Region Transfer, Ocean Desalinization, and 
Antarctic Ice Transport and eventual assignment of financial 
and economic resources to such projects. Inter-region 
transfer involving the movement of water from assured, 
excess supplies, for example, in the Northwest Territories 
in Canada to the Calif ornia-Nevada-Arizona area in the 
United States, would involve immense cost, formidable 
engineering difficulties and long time periods. Ocean 
desalinization, especially with the use of nuclear or solar 
technology, offers substantial opportunities for support 
of economic activity although at immense expense, with 
potential physical dangers and extended time of development. 
The transport of Antarctic ice and the melting of it in 
arid and semiarid coastal regions, although a relatively 
new idea and largely untried, offers possibilities for the 
support of economic activity. The transport of Antarctic 
ice presently appears to be comparatively inexpensive, 
involves the use of existing technology, and is capable of 
implementation in the near future for areas south and possibly 
even north of the equator. 

In arid regions, e.g., Sahara, Saudi, Atacama, Peruvian, 
Kalahari, and Australian Deserts, most parts of which lie 
far from adequate rainfall areas and desert rivers, the 
expansion of civilization will be dependent ultimately on 



ocean desalinization and the transfer of Antarctic ice. In 
semiarid regions, e.g.. South Africa, angola, Argentina, 

Iran, Pakistan, California, Mexico, Australia, etc., economic 
expansion and the support of larger populations will depend 
on exploitation of the local rainfall in combination with 
ocean desalinization, the transfer of Antarctic ice, and the 
expanded use of adjacent but distant stored rainfall. The 
western half of the United States, except for the relatively 
small, rain-rich Pacific Northwest stands as a monument to 
agricultural opportunity which is dependent upon the availa- 
bility of water. The imposing total of 15 western states 
is predominately semiarid or arid and dependent upon irri- 
gation works and scanty rainfall for the sometimes intense 
but mostly modest or nonexistent development of agriculture.' 1 ' 
Additional immense regions of the globe in similar favorable 
latitudes lie agriculturally dormant because of the water 
deficit and high evaporation and transpiration character- 
istics of the areas. 

In arid and semiarid areas lying south of the equator 
and within reasonable proximity to the oceans, e.g. , parts of 
the Republic of South Africa, scientific theory and experi- 
ment support a view that the transfer of Antarctic ice is 
limited more by the problems of potential resettlement and 



^Arizona, California, Colorado, Idaho, Montana, Nevada, 
New Mexico, North Dakota, Oklahoma, Oregon, South Dakota, 
Texas, Utah, Washington, Wyoming. 



the cost of irrigation works than by limits on the ice 
transfer technology. Given the relatively short distances 
from the tabular icebergs, and the comparatively cool 
intervening water and air temperatures for movement, the 
NPS scientists feel that scientific feasibility is close to 
being established for movement from Antarctica to South 
Africa, Argentina, Atacama, and Australia. Opportunities 
exist for the United States Government to increase its 
influence peacefully through the use of superior technology 
in Southern Africa, South America, and Australia. Engineering, 
economic, and political studies will be required to guide 
U.S. government policy makers in the exploitation of the 
scientific opportunities for the settlement of populations 
and the .expansion of agriculture and industry in previously 
excluded arid regions. 

The stages of transfer and use of Antarctic ice run in 
the following generalized pattern: 

Antarctic Ice : 

Stages of Transfer and Use 



South -• North 

Locate Breakout Transport Transport Terminal (Extract) 

\ \ / \ / \ / 'a 

Select Hookup Insulate Terminal (Dock) Distribution 



The stages through which an Antarctic ice transfer 
project must necessarily proceed are common for ice delivered 
to any location on the globe. The individual conditions of 



ocean routing, terminal docking, etc., may vary enormously, 
for example, between ice delivered to Southern California 
and South Africa, but the general process, or stages, of 
transfer will be virtually identical. With two noteworthy 
exceptions, the technology for each of the stages is largely 
in hand. Space photographic technologies are available to 
locate and select tabular bergs which have the sizes and 
shapes to meet other delivery criteria. The breakout of 
large glacial tabular bergs may not be necessary depending 
upon the time of the year, e.g., Antarctic summer, and the 
location of the glaciers relative to open water. The hookup 
and towing of the tabular bergs can be forecast to be 
engineering challenges of the first magnitude but surmountable 
with present technology. Ocean routing and terminal environ- 
mental impact should be able to be mastered with existing 
technology . 

The most imposing scientific challenge is the one of the 
melting of ice moving in seawater, and little theoretical or 
experimental knowledge is available on the subject. If 
rates of melt can be described adequately and predicted, the 
scientists can present the ultimate reality of ice transfer 
projects in terms of (1) the amount of melting and (2) the 
associated necessity of lack thereof for insulation. The 
most imposing potential problem in the transfer of Antarctic 
ice can be forecast to be the one of placing insulation on 
the ice and maintaining it there during the transportation 
to the terminal area. If the experimentally supported 
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theory shows that insulation is not necessary for the move- 
ment of tabular bergs to areas like Southern Australia, 
Argentina, and Southwestern Africa, the potential costs of 
movement would fall steeply, and the economics of transfer 
would be more attractive. The absence of insulation reduces 
costs and engineering difficulties to a point where the NPS 
scientists tentatively conclude that the transfer of ice to 
Australia, Argentina, and South Africa is technically and 
economically within the grasp of governments and the 
larger private businesses associated with those areas. The 
perception of the necessity for agricultural expansion on 
the parts of the governments and the economic atractiveness 
on the part of private business remain to complete a picture 
of Antarctic ice capable of being transferred to southern 
regions in the immediate (one to three years) future. 

The transfer of Antarctic ice to regions north of the 
equator present greater cost and engineering problems largely 
because of the problems of melting and the prima facie 
necessity for insulation. No realistic schemes presently 
exist for the insulation of the large tabular bergs in terms 
of physical or chemical barriers between ice and seawater. 

The NPS scientists, based on recognition of the fact that 
increased water temperature exponentially increases the rate 
of melt of ice, are gathering experimental data both in 
Monterey Bay and in the NPS mechanical engineering laboratories 
on the cooling effects of ice moving through water. A 
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"plume" of cool water is assumed to exist behind ice in 
water which could reduce water temperatures to the degree 
that other ice bathed in the plume would experience signi- 
ficantly reduced melting. The potential insulation scheme, 
in effect, would involve either (1) towing several bergs 
in column with a small "sacrificial" berg in front or 
(2) towing several bergs in column while removing the lightly 
compacted snow and ice from the top of the leading berg and 
depositing it in front of the leading berg as a coolant for 
the most exposed ice. 

Several challenging but surmountable problems exist at 
the terminal area for the tabular bergs. The ice upon 
arrival will have imposing drafts of approximately 200 - 500 
feet depending upon the distance and velocity of the tow 
and the water temperatures encountered. The immense drafts 
of the bergs will demand carefully selected terminal areas 
in order to place the ice conveniently close to the receiving 
water system and reduce the impact of adverse sea conditions 
on the floating reservoirs of fresh water. Ice, for example, 
delivered to Southern California could be "docked" immediately 
south of Santa Cruz Island, 20 miles from the California 
coast protected from the prevailing northwesterly winds and 
representing a reduced impact on the local environment. The 
most difficult engineering problem at the terminal location 
would be the one of establishing a barrier between the tabular 
bergs and the seawater in order to extract the most fresh 
water possible from the ice. Given the static situation of 
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the bergs at the terminal location, one can envision 
engineering solutions ranging from the production and trans- 
fer ashore of ice water slurries through the static wrapping 
of the ice in a waterproof plastic barrier and the extraction 
of the fresh water melt. 

If the technology can be developed to transfer Antarctic 
ice at reasonable cost to the coasts of arid and semiarid 
regions, immense opportunities would accrue to governments 
with the resources to exploit the technology. A government 
like that of the United States could work in a cooperative, 
positive way with the government of the Republic of South 
Africa to open up approximately 50,000 square miles of 
territory to agriculture and grazing in the more arid sections 
of the Cape of Good Hope and Orange Free State Provinces. 

The projects would be based ultimately on the technology 
necessary (1) to extract the maximum water from the existing 
rainfall and (2) to transfer and distribute additional water 
from Antarctic tabular bergs. The financial support for the 
water supply and agricultural project would involve invest- 
ment in a water distribution system and encouragement for 
settlement in the new territories. Such financing might be 
accomplished through a combination of U.S., Western European, 
South African government, World Bank, and private capital. 

The governments involved in cooperating in such a venture 
would be presented with opportunities to effect change in the 
hypothetical case of South Africa, without recourse to the 
hackneyed, aging and brutal techniques of insurgent and 
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counter-insurgent warfare, without the advantages accruing 
in insurgency situations to the weaker and less developed 
Eurasian communist countries, and with the strong possibility 
of effective change and political accommodation based on 
rising standards of living within the Republic of South Africa. 
The United States, in particular, with its advanced technology, 
large financial resources, powerful and advanced allies, and 
traditional predominance at sea is faced with the opportunity 
to use its natural strengths to support successful political 
accommodation and settlement in Southern Africa. 

In other parts of the world, in situations less traumatic 
than in South Africa, the U.S. has the opportunity through 
the exploitation of its technology to intensify agricultural 
activity in some areas and open up agricultural and asso- 
ciated economic activity in others. On the west coast of 
North America, the California and U.S. Governments are faced 
with a three-year drought crisis (1975, 1976, 1977). The 
crisis may require the positive supply of water in addition 
to the present negative provision of relief funds in order 
to master the situation and continue to support the existing 
overpopulation of California relative to scarce and ultimately 
limited water supplies. In the absence of the timely availa- 
bility of water from inter-region transfer and nuclear 
deslinization , Antarctic ice transfer probably will emerge 
as an imaginative, relatively low-cost alternative for the 
provision of fresh water to areas south of the equator. 
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With knowledge of the processes of melting, scientists, 
engineers, and governments could provide water to arid and 
semiarid regions north of the equator. 
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APPENDIX B 



SOLUTION OF TRANSIENT THERMAL CONDUCTION EQUATION 

The comparison between the theoretical curves and 
experimental data points contained in Chapter VI, Figures 
100 - 105 were obtained utilizing the computer program, 

"Reduce Data on Temperatures Internal to Ice", Figure B-l. 

The terminology used for the text and computer programs is 
listed in Table B-I. 

The computer program, "Select a best S", Figure B-2, 
was utilized to assist in obtaining values for thermal 
diffusivity, ice regression rate, and X Q for program Figure 
B-l. The value "S" is found by subtracting the experimentally 
determined values from the calculated values, squaring the 
difference and summing over the set of experimental data 
points. The theoretical function, which provides calculated 
values, was varied by changing the parameters a, X q , and V g . 

A minimum value of "S" represented the best set of above 
values . 

After experimenting with the aforementioned parameters, 
a minimum value of "S" of 1329 was attained using V g = 2.2 x 10 
and a = 3.0 x 10 these values of V g and a are utilized 

in the program shown in Figure B-l. 

The theory used to generate computer program 1 was the 
solution of the heat conduction equation as follows: 
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TABLE B-I 



NOMENCLATURE 



Symbol in 
Computer Program 
HP 9831 


Symbol 

in 

Text 


Definition 


Units 


H 


H 


heat of fusion for ice 


joule/kg 


Q 


q 


heat flux to melt ice 


joule/m^ sec 


A/Al 


a 


thermal diffusivity 
of ice 


2 . 

m /sec 


K 


K 


thermal conductivity 
of ice 


joule/msec 0 


V 


V s 


ice regression rate 


m/sec 


XO 


X o 


distance from edge 
of ice to first 
thermocouple 


meters 


Tl 


T 

s 


surface temperature 
of ice 


°C 


T0 


T. 

1 


interior temperature 
of ice 


°C 


X9 


AX 


thermocouple spacing 


meters 


I 


I 


time intervals (1-6) 


minutes 


P 


P 


density of ice 


kg/m^ 


W 


W 


V s /a 


1/meters 



99 



SEAWATER 






'W 



ICE 



/ 



■>X 



x=o 



/ 



x=x 



An energy balance at the water-ice interface from above is : 



q = Q + q . = 

^w w 



HpV s ' <k 3X>X-X (B - 1) 



Of the heat transferred from the water the amount of heat 
which causes melting is: 



Q 



H p 




HP V S 



(B-2 ) 



For the purpose of analysis, the temperatures within the 
ice are modeled using the one-dimensional unsteady heat 
transfer equation as follows: 




1 3T 

a at 



0 



(B-3) 
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The following assumptions were made: 



1) q = constant 

w 

2) Steady state melting occurs 

3) V s = constant = velocity of surface of ice 

The boundary conditions at the water-ice interface are as 
follows : 



which mathematically is as X °°. 

Define a variable which involves both space and time 



Note that the solution of eq. (B-3) will use £ as an 
independent variable. Consequently the temperature profile 
propagates with the same shape into the ice as the ice 
melts. Combining equations, one can derive the following: 




X = X 




(B-4) 



x 



T [X s ,t] 



T 



s 



(B-5 ) 



The boundary condition for the ice interior is: 



T [°°, t] 




(B-6 ) 





(B-7) 



-iV T i lf " - ^'WV 1 



0 



(B-8 ) 
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where the function "f" is given by: 



T(X,t) = T - [T -T. ] f(£) 

o SI 



(B-9 ) 



We now define: 



ip = f ' 



(B-10 ) 



V V 

f" + — f ’ = ib’ + — ih 

a a 



= 0 



gf + — ♦ = 0 

d£ a 



V 



d \jj _ _s 






dC 



a 



ip = ae 



-V s 5/a 



Taking the derivative of f with respect to E, and then 
integrating to solve for f as follows: 



df - v s i/a 

dT - ae 



, , - v s 

f = J ae 



dC = - a— 



v s y v ^ /a + h 

- a — e + b 

a 



c = - a V / 



a 



-V_ C/a 

f = ce + b 



(B-li: 
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Utilizing boundary condition (6) and combining eq. (B-9) 
and (B-ll ) we obtain the following: 



T[»,t] = T - [T -T. ] [ce"°° + b] = T, (B-12) 

O O -I- 

b = 1.0 

Utilizing boundary condition (5) when X = X g we obtain the 
following: 



C = I = X - V t = 0 

s s s 



(B-13) 



T [X ,t] = T c = T - [T -T • ] [ce"°+l] 

o ^ o o o 



c = -1 



It therefore follows that: 



TCX^t] = T s - CT S -T i ] [1 - e 



“V £/ ot 



] 



(B-14 ) 



C = x - V t > 0 

s — 



Utilizing boundary condition (4) we obtain the following: 



111 

3X X=X 



fV T i 1 [ -V“ ] e - V /a at 



_ 



V 

- — [T -T. ] 
a s l 



(B-15 ) 
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We note that when X-X , £ = 0, therefore: 

-k V • 

[T -T . ] = HpV - q 

a si K s 

q w = Hp v s + k/a V s [T s -T i ] 

The final result is as follows: 

[T -T-] 

q w = V s [Hp + k - - - -- -- - ] (B-16 ) 
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10 PRINT "REDUCE DATA ON TEMPS INTERNAL Tu ICE 
20 PRINT "TRACK 0« FILE 1" 

36 REM k > sub s is u! alpha is a. 

4u REM % sub s is tl! t sub ice is t0. 

50 DISP "PLOT ?YE3=1 ? NO=0" 5 
60 INPUT GO 
70 DISP "INPUT V"; 

80 INPUT V 
90 DISP "INPUT Al" i 
100 INPUT hi 
110 M 0=720 

120 DISP "INPUT XO"! 

130 INPUT X0 

140 DISP "INPUT I"; 

150 INPUT I 
160 W = V / A 1 
170 MS =300 
ISO Tl =0 
190 T 0=-2 1 . 7 

200 PRINT " SUPF TEMP IS "Tl" DEG C. ICE INTERIOR TEMP IS "TO" HE 
210 PRINT " C." 

220 PRINT "V="V" THERMAL D I F F U S I V I T V " A 1 " M T 2 ' S EC. " 

230 PRINT "X0="X6 
240 REM g reek 1 is z. 

2 5 0 P E M s 3 i n i v t e t e r s . 

260 3 3= X 9=0. 00158 

270 DISP "DATA TABLE n VES=1 * NO=0"S 

230 INPUT Z6 

290 PRINT 

30O FORMAT 4F16.4 

310 IF II THEN 360 

320 A =70 

330 T8=0 

340 PRINT "TIME T8="T3 
350 GOTO 390 

3 60 A = X 0 - V * '• M 0 + CI-D - M 9 ) 

370 T8=M0+'. I-1^M9 

330 PRINT "TIME T8="T3. 60 

390 B=A+4+X9 

400 IF 00=0 THEN 490 

410 DISP "INPUT V AXIS SCALE C"3 

420 INPUT C 

430 DISP "INPUT V AXIS SCALE D"J 

440 INPUT D 

450 P=B-A 

460 PRINT " R= ""R 

470 N = <: B-A * -- S3 

430 PRINT " N= " N 

490 PRINT " TEMP C TEMP F XI* METER" 

500 IF 00=0 THEN 540 
510 SCALE A ♦ B < C * D 
520 X AX I S D ? < B-A 1 4 
530 VAX IS H» •• D-C ■ 5 

549 N = *:.B-A>/S3 

550 R=B-A 

560 FOP Z=A TO A+R STEP' RX N 
570 T=Tl-'-Tl-T0>^‘ 1-E::P(.-W-2 > • 

530 T3=1.8*T+32 
590 IF Z6=0 THEN 620 
600 WRITE ( 2 • 300 '- 1 T * T3* Z 
610 IF 00=0 THEN 630 
620 PLOT Z,T 
630 NEXT Z 

640 IF 00=0 THEN 740 
650 FOR J=1 TO 5 
660 08= D + D+C )■■ 5 
670 ''AXIS 08 
630 NEXT J 
6 9 0 F 0 R r = 1 T 0 4 
700 07 =A+K ^<B-A;< 4 
710 VAX IS 07 
720 NEXT K 
730 PEN 
740 STOP 
750 END 



B-l. Ice interior temperature program 
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tor pro ? raft 1 



10 FP I NT " T h i i p r o j r a ft w ill s e 1 e c t a b e £ t. S 
20 PRINT "Pro 3 raft stored in track o< tile 2." 
30 DIN ::C5 ]jNC5] 

4Q DIM YC5i51«TC5»5>Z[5s5] 

50 GOTO 530 
bU M 9=300 

7 0 X 9 - 1 •••■ '• 39. 37^16 ) 

80 DISP "INPUT R"; 

90 INPUT H 

100 DISP "INPUT .:o m ; 

110 INPUT X0 

120 DISP "INPUT V"; 

130 INPUT V 

140 DISP "TABLE WANTED YES=1 NO=0' ; 

150 INPUT US 

160 F O R M H T 2F4 . U < F7 . 3 1 F9 . 0 ' j.F'9 . 1 < F 7 . F 3 . U 
170 IF QQ=0 THEN 190 

ISO PRINT "N JiX'l JmM'- I' ' T « > J *,Y'* I, J ■ ,$" 

190 DISP "INPUT M0 " i 
200 INPUT M0 
210 T 1 = 0 
220 TO =-10 
230 W=V/h 



240 S=0 

250 FOP J=1 TO 5 
26u : :c J > ’0+' J-i •*, ;s 
270 FOR 1=1 TO 5 
280 MCI>M0+I-M9 

29u zc i ■« j 3=:;c j i-vwic i J 

300 tc 1 1 j ] = t i - 1 t i - t o ■ ~*. i-e::p* -w-:t i - j a 

310 81 = ': T C I, Jl-YC I ' J ] 1 1 2 

320 S=S+3 1 

330 IF Q8=0 THEN 350 

340 WRITE ■ 2 « 1 60 :■ I ' J - . C J1<MC I J i TC I i J 1 < TC I « ,• 3 * 2C I i 
350 NEXT I 
360 NEXT J 



370 


PRINT "SUM S IS" 1 


s 


380 


PRINT 




390 


PRINT 


310 = " 1 


400 


FP I NT 




410 


S = u 




420 


GOTO 80 




430 


PRINT “ 


I 


440 


FOP 1=1 TO 5 




450 


FOR J=1 TO 5 




4 6 O 


WRITE 1 2 ’ 500 * I < J 


i YC I- 


470 


NEXT J 




430 


PRINT 




490 


NEXT I 




500 


FORMAT 2F 12. Oi FI 


S • o 


510 


STOP 




520 


END 




530 


YC 1 i 1 >*.-2.5-32.) 


1 . 3 


540 


YC 1 1 2 ]=*.:- 3- 32. 1 1 




550 


TC li 3 ]=*: - 3. 5- 32 


1.3 


560 


YC 1 i 4 ]=* -3.5-32 


1.3 


570 


'TC 1 1 5 1=' -3.5-32) 


1.3 


530 


'TC 2i 1 ]=< -0.7-32') 


1.3 


5 9 0 


Y'C 2 1 2 ] = * -1-32 1 1 


. 3 


6O0 


YC 2 j 3 ]= * -1.4-32) 


1 . 3 


6 1 O 


YC 2« 4 >c -1.5-32) 


• 1.3 


620 


YC2> 5X-1.6- 32) 


. 1.3 


630 


'TC 3 1 1 >*.2.3-32)/ 


1.3 


6 4 0 


YC 3i 2 3=*' 1.9-32)/ 


1.3 



650 YC 3-3]=* 1.5-32')- 1.3 
660 YC 3 ’4 ]= '. 1. 3- 32. -. '1.3 
670 YC 3? 5 3=' 1. 1-32 - 1.3 
630 i'C 4 . 1 >.. 8,4-32..'/ 1.3 
6 9 O Y C 4i 2 ]=<7. 7-32 1 1.8 
7O0 YC 4 1 3 ]= ' 6 . 3- 32 ' 1 . 3 

710 YC4,4]=> 6.5-32 1 1.3 
720 YC 4 1 5 ]= 1 6-32 1 1.3 
730 YC 5 ' 1 ]=*-22. 5- 32 ; 1 1.8 
7 40 YC 5 1 2 1= < 21 -32 - 1.3 
750 'rC 5' 3 ]=■. 19. 3-32 ■ 1.8 
760 YC 5* 4 ] = < 13. 7-32). 'I . 3 
7 7 o YC 5< 5 ]=< 17. S - 3 2 ) "1.3 
780 GOTO 60 



M 0 = " M0 



J ]< S 



B-2 . 



Program for minimum 1 S ' 



value . 



106 











REYNOLDS 
NUMBER 
- (LOG 10) 




11 


- 






10 


“ANTARCTIC 

□ 






9 










FL0W : TURBULENT 




FLOW 5 TURBULENT 


8 


- BOW WAVE 5 

NOT ATTACHED 




BOW WAVE 5 

ATTACHED 


7 


- NPS Re *" 2 
BAY a 




Region 3 


6 








5 


FLOW 5 LAMINAR 
“NPS BOW WAVE 5 
LAB ONOT ATTACHED 




FLOW 5 LAMINAR 
BOW WAVE 5 

ATTACHED 


4 


Region 1 




Region 4 


* 


1 L_ 




_J I 



0 0.1 0.2 0.3 0.4 



FROUDE NUMBER 



Regimes of flow and wave producing conditions as 
functions of Reynolds and Proude numbers. 

1. Regimes of flow and wave producing conditions as 
functions of Reynolds and Froude Numbers. 
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NPS water tunnel 
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109 



NPS mobile water table. 
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NPS water channel 
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NPS water channel. 
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Progressive regression of second laboratory model 
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Progressive regression of streamlined laboratory model. 




8. Effectiveness of heat transfer. 
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Flow reattachment with accompanying large heat flux. 
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10. Trapped vortex filament. 
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11. Ice model freezing container. 
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Top view of ice model freezing container 



GLACIAL ICE TOWING BOLLARD 



36" 



-o- 



I o 
!o 



J 



STANDARD 2X4 
18 STANDARD 4X4 H 

(OAK WOOD) 

X 2^2 lag screw 
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2 
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^ 7 II 

*3 ► 



o 

o 
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-<L 



^ 74 " plywood 



if* 1 ** @ 



FRONT VIEW 

Note’- Treat assembly with 
Thompson water seal 
or equivalent 

^2X2 metal angle 
bracket 



^gcarriage bolt 
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14. Premature disintegration of ice model. 





o 
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GLACIAL ICE CORE BOX 




16. Core box installation. 
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17. Ice model underway with core box installed. 
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18. Ice model with towing bollard integral with core box. 




19. Scale used to weigh ice and container. 
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20. Positioning of ice container near water and towing vessel. 
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21. Ice container being lowered into the water. 
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22. Seawater flooding into ice container. 
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23. Ice model separation from container. 
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24. Ice model underway. 
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Coast Guard vessel with model in tow 
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26. Boston whaler with model in tow. 
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27. Tow rope and instrumentation wire arrangement. 
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meter and readout assembly 
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29. Ice model positioned for recovery. 
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30. Recovered ice model clear of water 
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31. Ice model in pier; Test 
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32. Ice model underway; Test 
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33. Premature disintegration of ice model; Test 
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34. Ice model underway; Test 






m 
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35. Ice model recovery; Test 
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model recovered 
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37. Ice model on pier; Test 
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38. Ice model bow erosion; Test 
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39. Ripples along the complete port side of the ice model; Test 
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40. Detailed view of port side, rear of ice model; Test 





GLACIAL ICE THERMOCOUPLE POSITIONS 
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41b. Thermocouple with digital readout. 
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arrangement; Test 




I 
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Ice model underway; Test 
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Forward section of ice model underway; Test 
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Bow undercutting of ice model; Test 
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46. Ice model recovery; Test 
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Recovered ice model with exposed thermocouple 
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Ice model on pier; Test 
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Detailed view of ice model layering and ripples; Test 
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50. Ripple measurement technique; Test 
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51. Ice model prior to launch; Test 



GLACIAL ICE THERMOCOUPLE POSITIONS 
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Forward ice 
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52. Thermocouple arrangement; Test 5. 
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160 



Cross section view of pontoon bridge. 
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54. Thermocouple installation on pontoon bridge; Test 
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55. Ice model and pontoon bridge underway; Test 
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Ice model and pontoon bridge towing arrangement; 
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57. Ice model on pier; Test 
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Rear view of ice model ripples; Test 
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59. Rear view of ice model layering; Test 
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60,. Thermocouple arrangement; Test 6. 
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61. Ice model prior to launch; Test 
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Forward section separation; Test 







I 
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63. Forward section separation; Test 
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64. Rear view pushing arrangement; Test 
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Frontal view pushing arrangement; Test 
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66. Forward section of ice model; Test 
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67. Ice model positioning for recovery; Test 
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68. Frontal view of recovered ice model; Test 
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69. Reference grid on ice model for ripple measurement; Test 
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70. Linear motion position transducer used for ripple 
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71. Strip chart recording of ripples; Test 
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72. Ice model immediately after insertion in water channel. 







180 



73. Ice model after two minutes in water channel showing dye injection. 
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74. Tandem models in water channel immediately after insertion - separation 
distance = W/3. 







in 
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Tandem models after one minute in channel at W/3 separation distance 




183 



76. Tandem models after two minutes in channel at W/3 separation distance. 
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77. Forward ice model bottom view. 
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Rear ice model bottom view. 
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79. Dye accentuated flow visualization around tandem models at 1/3W 
separation distance after two minutes in channels. 
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80. Hookup of tandum ice models in Monterey Bay. 
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81 . Side view of tandem ice models underway. 
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82. End view of tandem models underway. 
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83. End view of tandem models underway oscillation left. 
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84. End view of tandem models underway oscillation right (Note drag rope 
from rear model) . 
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85. Bow view of forward tandem ice model on pier. 
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86. Bow view of rear tandem ice model on pier. 
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LOCAL ICE ABLATION RATE, 




90. Local ice ablation rate: 2 October 1977. 
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CURVE PREPARED BY 0 M. GRIFFIN, NRL (1/27/78) 
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Turbulent flow Stanton number vs. Reynolds number 
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92. Turbulent flow Stanton number vs. Prandtl number. 
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BOUNDRRY LRYER PROFILE 
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93. Boundary layer profile, 6 November 1977. 
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RFTER BOUNDARY LAYER PROFILE 
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95. After boundary layer profile, 10 December 1977. 
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96. Ice interior temperature profile, 6 November 1977. 
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97. Histogram of ripples along ice, 3 October 1977. 
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Histogram of ripples along ice, 10 December 1977. 
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99. Histogram of ripples along ice, 27 April 1978. 
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INTERVAL 2 (TIME=17 MIN. 
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DISTANCE IN ICE (METERS) 
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